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ABSTRACT 

We report on multi-epoch optical and near-infrared spectroscopy around the first overtone ro- 
vibrational band of CO in the pulsating yellow hypergiant p Cas, one of the most massive stars 
in the Galaxy and a candidate SN II progenitor. We argue that the double cores of the CO absorption 
lines, that have previously been attributed to separate circumstellar shells expelled during its recurrent 
outbursts, result in fact from a superposition of a wide absorption line and a narrow central emission 
line. The CO line doubling returns over subsequent pulsation cycles, where the superposed line emis- 
sion assumes its largest intensity near phases of maximum light. We find that the morphology and 
behavior of the CO band closely resemble the remarkable "line-splitting phenomenon" also observed in 
optical low-excitation atomic lines. Based on radiative transport calculations we present a simplified 
model of the near-IR CO emission emerging from cooler atmospheric layers in the immediate vicinity 
of the photosphere. We speculate that the kinetic temperature minimum in our model results from a 
periodical pulsation-driven shock wave. We further discuss a number of alternative explanations for 
the origin of the ubiquitous emission line spectrum, possibly due to a quasi-chromosphere or a steady 
shock wave at the interface of a fast expanding wind and the ISM. We present a number of interesting 
spectroscopic similarities between p Cas and other types of cool variable supergiants such as the RV 
Tau and R CrB stars. We further propose a possibly common mechanism for the enigmatic outburst 
behavior of these luminous pulsating cool stars. 

Subject headings: infrared: stars — line: profiles — shock waves — stars: atmospheres — stars: 
oscillations — supergiants 



1. INTRODUCTION 

p Cas (HD224014) is a naked-eye variable supergiant 
star (V = 4 m .5, K = 2 m ) with a long history of photo- 
visual photometric and spectroscopic observations. It 
was observed with the first near -infrared (NIR, 1-3/xm) 
spect rographs in the late 1970s l|Lambert. Hinke fc Halll 
Il981|) , but has been somewhat neglected since then. This 
paper discusses recent spectroscopic monitoring of the 
NIR spectrum of this remarkable hypergiant star, reveal- 
ing a wealth of remarkable temporal spectral features 
whose interpretation remains elusive. For example the 
cores of numerous strong absorption lines always appear 
double. The apparent line core splitting results from a 
static central emission component. However in certain 
variability phases p Cas also reveals a prominent opti- 
cal and NIR emission line spectrum that can indicate a 
common physical origin with the central emission com- 
ponents of split absorption lines. An important question 
to further address the peculiar nature of p Cas' emission 
spectrum is therefore if the central emission lines and 
the prominent emission line spectrum form in similar cir- 
cumstellar conditions (e.g. at similar distances from the 
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surface of the hypergiant with comparable geometric den- 
sity distributions). Detailed spectroscopic studies of the 
circumstellar environments of cool hypergiant stars are 
scarce, although they are absolutely required to prop- 
erly understand the physical mechanisms that can ex- 
plain their uncommonly large observed mass-loss rates. 
For example, can the large mass-loss rates observed in 
F- and G-type hypergiants, such as p Cas (that are too 
cool to accelerate the radiatively line-driven winds of lu- 
minous hot stars), result from atmospheric pulsations, 
since their extended winds are too warm to be accel- 
erated by dust-driven mechanisms commonly proposed 
in red supergiants? Can the observed properties of p 
Cas' variable emission line spectrum be reconciled with 
the pulsation kinematics of the extended atmosphere and 
with its uncommonly large mass-loss rates? 

With a spectral type of F8-G2 IaO (T cff = 6500- 
7200 K), and a luminosity of log L*/Lq ~ 5.7 (i?* 
~ 400-500 R©), p Cas is one of a small sample of 
known yellow hypergiants in the Galaxy. In quies- 
cent variability phases the star looses mass at a rate 
of few 1O _4 M0, which is intermediate for cool hyper- 
giant stars, but can increase by two orders of magni- 
tude during rare outburst episodes llLobel et alJ 119981 
120031). Evolutionary tracks (e.g. iLeieune fc Schaererl 
(2001)) a nd possible memb ership to the Cas OB5 as- 
sociation (Humohrcvs 1978) indicate a stellar mass of 
~40 Mq and possible ages of 4-6 Myr. Observations of 
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CNO-proces sed abundances in its extended Na-enriched 
atmosphere ijTakeda fc Takada-Hidailll994D signal that 
this yellow supergiant has gone through dredge-up in the 
red supergiant stage, and is presently evolving bluewards. 
The star is about to cross the so-called atmospheric in- 
stability region in th e upper HR-diagram: "the Yellow 
Evolutionary Void" l|Nieuwenhuiizen fe de Jagerl fl995'). 
where due to the enormous atmospheric scale height of 
these massive evolved stars, coupled with a large tem- 
perature in the shell-burning layers, the effective surface 
acceleration essentially vanishes, producing dynamically 
unstable stellar atmospheres. 

p Cas's atmospheric instability is manifested in both its 
quasi-periodic photometric behavior and unusual spec- 
tral variability. It normally pulsates semi-regularly (Srd) 
with P = 320-500 d, and an amplitude of ± m 2 
(V band) (IZsoldos fc PercvllT99ll lArellano FerrollT^ . 
Once every 20 to 50 years, however, this hypergiant goes 
into an outburst by dimming more than a full visual mag- 
nitude, decreasing the T e g by more than 3000 K, and by 
ejecting a massive gas shell. The most recent outburst 
event occurred in 2000-01 which was carefull y monitored 
using high-resolution optical spectroscopy l)Lobel et al.1 
2003). In this paper we present the NIR (2 pm) spec- 
trum from that exceptional epoch in H3.ll 

The spectrum of p Cas outside outburst is pecu- 
liar. Optical studies reveal that the spectral lines are 
very broad due to large atmospheric micro- and macro- 
turbulence velocities (11 and 21 kms -1 , respectively). 
The detailed profiles of certain lines are also rather com- 
plicated. Besides the single-component weak absorption 
(" normal" ) lines that are entirely formed in the super- 
giant's lower photosphere, the low-excitation (xiow ^ 3 
eV) absorption lines are either cyclically or permanently 
split (depending on the line oscillator strength within a 
multiplet), having one component always blue- and the 
other one red-shifted relative to the stellar radial veloc- 
ity of — 47±1 kms -1 . A number of low-excitation lines 
appear in emission above flux level of the local contin- 
uum around variability phases of maximum light (e.g. 
when T e g becomes largest), while others, like the for- 
bidden [Call] lines, peak always prominently above the 
local continuum level throughout a complete variability 
cycle. In addition, many medium-strong and strong ab- 
sorption lines develop far violet extended wings during 
certain phases of fast wind expansion. They can some- 
times also display a rather "bumpy structure" in the far 
short- wavelength wing profile of split absorption lines. 

Despite the fact that p Cas's peculiar spectrum was 
immediately recognized almost a century ago, its pre- 
cise physical origin remains heavily debated. This is no 
surprise given that it requires an appreciable amount of 
observing time to acquire the substantial body of simul- 
taneous spectroscopic and photometric data necessary to 
reveal a correlated pattern of behavior. The earliest in- 
terpretations of the peculiar line splitting phenomenon 
attribute it to two detached circ umstellar shells previ- 
ously ejected by the superg iant ( Bidclman & McKcllar 
119571 lSarge"rr3ll961llG~esickilll992|) . Various authors re- 
marked that such shells possibly comprise two atmo- 
spheric layers flanking an outwardly propagating shock 
wave, a mechanism currently used to explain the ab- 
sorption line splitting observed in a number of less- 
luminous variables such as RR Lyr, W Vir, RV Tau, 



and Miras (e. g.. IChadid fc GilleU Il99?l IWhltnevI 
T962t ICillet et a.UI1 989HHmkle. Hall, fc RidgwayHl982ir 
Shcffcr & La mbert! (|1986f> added support to the shock 
wave hypothesis for p Cas by showing that the line core 
doubling correlates with the pulsation phase, because 
the splitting is most prominent around maximum light 
when the photosphere begins to expand. They also ar- 
gued that the line core doubling first starts in the high- 
excitation lines, as if a shock wave were propagating out- 
wards crossing vari ous atmospheri c regions with decreas- 
ing xiow However, Shcffcr (1993) also considered a dif- 
ferent explanation for the split lines of p Cas. Instead of 
two separate absorption components, the splitting could 
be interpreted as a superposition of a single broad ab- 
sorption line and a static central emission line. This 
proposition was based on the observation of the cores 
of some split lines that become so intense at times that 
they peak above the local level of the continuum, some- 
times without any accompanying absorption, which also 
rules out a P Cygni-type line formation region in p Cas. 
These effects were also noticed for the slightly warmer 
hypergiant HR 8752. The emission can be better de- 
tected aga inst the weaker absorptio n spectrum of this hy- 
pergia nt l)Sheffer fc Lambert! Il987j) . iSheffer fc Lambert! 
l)1987|) hence proposed that the origin of the variable 
emission lines could be in a putative extended chromo- 
sphere around these very luminous cool stars. 

Further steps at discarding the "shell hypothesis" for 
the form ation of split absorp tion lines in p Cas were pro- 
vided in lLobel etaD l|1998f) . They showed that 1) the 
observation of permanently split absorption lines in the 
optical spectrum, with components remaining at a con- 
stant radial velocity, would require two absorbing 'shells' 
moving with positive and negative radial velocity dur- 
ing the complete variability cycle. Two separate shells 
permanently propagating in opposite directions is physi- 
cally excluded; 2) the hig h-excitation lines considered in 
ISheffer fc Lambertl <)1986j) do not necessarily form deeper 
than low-excitation lines in the extended atmosphere of 
p Cas, and therefore do not provide additional 'evidence' 
for an outwardly propagating atmospheric shock wave; 
3) the variations of the red- and blue-shifted components 
in the split absorption lines are not independent of each 
other. One line component always grows stronger while 
the other line component weakens. This effect can read- 
ily be mimicked by assuming a central static emission 
line superimposed on a stronger and broader absorption 
line that exhibits small Doppler displacements due to 
the photospheric oscillations. Further direct evidence 
against substantial circumstellar material around p Cas 
includes the non-detection by HST-WFPC2 of any sig- 
nificant structure at 0.1" or 310 AU » 150 R+ away 
from the hypergiant ijSchuster. Humphreys fc Marengol 
2006). Nevertheless, tenuous and temporal circumstel- 
lar material m ay still be present as ind icated by a weak 
mid-IR excess ijJura fc Kleinmannll990|) and the spectro- 
scopically observed mass ejections during recurrent out- 
bu rsts of the s t ar (e. g. optical TiO band observations 
bv lLobel et alJ l)2003f) and references therein to the two 
earlier recorded outbursts). 

In this paper we address the important question of the 
origin and nature of the line emission region in p Cas. 
In particular, what are its thermal conditions and what 
is the geometrical extension and location of this region? 
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Can the emission region and its line excitation mecha- 
nism be linked to the pulsations of the lower photosphere 
of the hypergiant? We discuss several important clues 
from NIR spectra of p Cas, but will also scrutinize new 
high-resolution optical and UV observations (of which a 
more detailed analysis will be presented elsewhere). 

While optical TiO bands have only been recorded dur- 
ing rare outbursts, the abundant CO molecule is how- 
ever always observed in the NIR spectrum of p Cas. The 
NIR spectrum was extensively monitored by Sheffer et 
al. (1993) with the FTS spectrograph on the Mayall 4m 
telescope (Kitt Peak) from 1979 to 1989. Their intense 
monitoring program revealed substantial NIR spectra l 
changes on a time scale of only months (Sheffer 1993). 
Although the "shell hypothesis" can be abandoned for 
the optical split absorption lines, these authors adopted 
it to interpret the CO line variability in the NIR spec- 
trum. Unlike the optical atomic split lines, the NIR CO 
lines require kinetic gas temperatures much smaller than 
those adopted in classical photospheric models for F- and 
G-type supergiants (fo r a comprehensive CO Grotria n 
diagram see Fig. 1 of Bie eing. Rieke. fc Riekd l|200^ ). 
The low kinetic gas temperatures (< 2000 K) required 
for CO molecular formation at these early spectral types 
either require a temperature reversal in the upper at- 
mosphere (a model we explore in i)4.2f> . or must assume 
the presence of a much cooler circumstella r gas envelope . 
The latter interpretation was adopted bv lShefferl lli"99"l 
by attributing the short- wavelength " shell" component of 
the split CO absorption lines to a moderate outburst of 
p Cas in 1986. They further assumed the detached shell 
interpretation to compute the total mass of the expelled 
gas layer, its turbulent velocity, and the atmospheric de- 
celeration during their monitoring period. 

NIR CO bands have been observed and interpreted as 
direct tracers of ejected stellar materi al in a number of 
other F-G supergiants. For example, lOudmaiier et al.l 
ill 9951) discuss the m o lecula r winds of post-AGB stars, 
while iGonzalez et alJ (jl998) investigate dust formation 
close to the stellar surface of R CrB stars. Given the dif- 
ficulties to properly explain the properties of split atomic 
lines in p Cas, we argue that a "shell hypothesis" for CO 
line formation needs to be re-examined as well. The early 
investigations by Sheffer already indicated that despite 
the wide variety of morphologies exhibited by the CO line 
profiles, they can nonetheless be reduced to a narrow cen- 
tral emission line superimposed on a broad absorption 
line, yielding two adjacent (but therefore only "appar- 
ent" ) absorption line components. However, their study 
was not widely published or disc ussed, and only one FTS 
spectr um has been discussed bv lLambert. Hinke fc HalU 
IU981 . The major part of these exceptional NIR spec- 
trosco pic data were only presented as part of a Ph.D. 
thesis llShefferlll 993ft . 

We present a collection of NIR spectra of p Cas ob- 
tained with several telescopes in 2000-04. Our NIR moni- 
toring program is not as frequent as Shcffer's, but has the 
following benefits. First, we present spectral analyses of 
new NIR observations in combination with our ongoing 
optical hig h-resolution spectr oscopic monitoring of p Cas 
since 1993 l|Lobel et al.ll2003|) . Next, while Sheffer et al. 
studied NIR spectra of the moderate outburst of 1986, we 
investigate the major outburst event of 2000-01. Follow- 
ing the latter major outburst the long-term variations 



in Ha reveal the formation of an exceptio nal P Cygni 
type line profile since the spring of 2004 ifTohell 12004: 
lLobel et al.ll2005j) . Cool supergiant stars often exhibit 
secular variations on time scales of decades (for exam- 
ples see t!5.2.1J) and it is therefore important to continue 
accumulating spectroscopic observations of these unusual 
stars. 

Although our combined spectra span the entire range 
from optical wavelengths longward of Ha through the 
main NIR J, H, and K bands, up to L (4 (an), the 
present study primarily concentrates on the K band re- 
gion where the first overtone band of CO occurs (JJ. 
Our high-resolution CGS4 spectral observations of 2004 
September are motivated by the large variability of this 
molecular band observed in low-resolution spectra, and 
because it was detected in strong emission in 1998. We 
also draw attention to the L band region in the SpeX 
spectrum Ci'i.l\ because, to our knowledge, it is the 
first L band spectrum of medium spectral resolution pre- 
sented for any F-type supergiant like p Cas. 

The outline of this paper is as follows. First we de- 
scribe new optical to NIR spectra of p Cas obtained in 
1998-2004, covering roughly four pulsation cycles (SJ3J. 
Then we establish the presence of emission in NIR CO 
lines and investigate their profile changes with the stellar 
pulsations ([JSJ- We discuss a unified excitation mecha- 
nism that can produce emission in both molecular and 
atomic lines, considering several explanations based on 
detailed CO line modeling, combined with co-eval optical 
spectroscopy (2J. Finally, we compare our observations 
with the line splitting phenomenon reported in other cool 
star spectra. For example, we show that brightness de- 
clines of R CrB stars appear intriguingly similar to the 
outbursts of p Cas (fJHJ- The comprehensive approach 
of our NIR monitoring program enables us to investigate 
the remarkable evolution of the CO 2.3 pm band, and 
to provide credible evidence for the stellar pulsations as 
the dominant factor that determines the profile changes 
observed in NIR and optical low-excitation spectral lines 
(©■ 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. IR Spectroscopy 

We have collected a series of low (R ~ 300) to 
high (R ~ 40,000) spectral resolution NIR spectra of 
p Cas covering roughly four pulsation cycles in 1998- 
2004. Our observations utilize three instruments on three 
telescopes - the CorMASS echelle spectrograph on the 
1.9m Vatican Advan ced Technology Telescope in Arizona 
ll Wilson et al.ll200l|) . the SpeX ech elle spectrograph o n 
the 3m InfraRed Telescope Facility l|Ravner et alJ l2003). 
and the CGS4 spectro graph on the 3.8m Unit ed Kingdom 
Infra- Red Telescope ({Mountain et alJ fl990(l at Hawaii. 
Observations with the latter spectrograph were obtained 
in service mode. The datasets are complemented with 
two archival CGS4 spectra, one obtained in 1998, one 
pulsation cycle before, and the other one right in the 
middle of the grand outburst minimum of 2000-01. Table 
Q] summarizes all observing dates and the basic charac- 
teristics of the new NIR spectra. 

All NIR observations are performed in a standard fash- 
ion resulting in similar data reduction procedures. A 
series of successive A-B sequences (where A and B cor- 
respond to two positions of the star along the slit) are 
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obtained, both for p Cas and a bright telluric standard 
A-, F-, or G-type star. Exposure times range from 0.5 
to 2 s, usually restricted to the shortest time allowed for 
a particular instrument. The size of the entrance slit 
ranges from 0.3" to 2", the allowed minimum to attain 
maximum available spectral resolution. To suppress the 
brightness of p Cas before saturation a neutral density 
filter is used with CorMASS and for the low-resolution 
mode of CGS4 (both for p Cas and telluric standards). 

The flat-fielded A and B images are combined and mu- 
tually subtracted to perform dark and sky removal. The 
spectral orders are extracted, wavelength-calibrated, and 
divided by a standard star spectrum with removed intrin- 
sic hydrogen lines. Next they are trimmed and combined 
into a single NIR spectrum. The resulting spectrum is 
then multiplied by a black-body spectrum corresponding 
to the spectral type of the standard. The wavelength 
calibration of the CorMASS spectra is performed using 
observations of the planetary nebula NGC 7027 obtained 
during the same night or the night before. The calibra- 
tion of the high-resolution CGS4 spectra and the short- 
wavelength mode of SpeX employ Ne-Ar-Xe arc lamp 
exposures obtained immediately after science target ob- 
servation. The low-resolution CGS4 observations and 
the SpeX long-wavelength mode calibration both em- 
ploy these arc lamp exposures as well as extraction of 
sky lines from the target images. An accurate absolute 
flux calibration has not been attempted for the obvious 
reasons of intentional flux suppression and strong seeing- 
induced variations between subsequent exposures of this 
very bright star. 

A number of important differences in our reduc- 
tion procedures using different instrument/grating set- 
tings are the following. The 40 1/mm CGS4 obser- 
vations require one additional step of correcting for 
a 2 to 3 pixel-wide wiggle in the flux level intro- 
duced by variations in the seeing conditions. It oc- 
curs when the detector is shifted in the dispersion direc- 
tion between subsequent exposures to improve the sam- 
pling of the resolution element on the detector array. 
The basic steps for reduction of the calibration frames 
and co-addition of the target frames are performed 
with the ORAC-DR pipeline for the CGS4 low-resolution 
spectra. The SpeXtool software package is used for 
the SpeX data which also performs the extraction of 
the final 1-d spectra l | Va cca 1 pushing fc Ravnerl 2003; 
ICushing. Va,cca fc RavneHl2004D . The CorMASS spec- 
tra are reduced with the standard iraf software pack- 
age imred.echelle. The SpeXtool routine Xtellcor 
provides a modified spectrum of Vega to perform tel- 
luric corrections in the SpeX spectra. For all other cases 
strong hydrogen lines in the standard spectra are inter- 
polated with a straight line and removed using the iraf 
task splot. 

2.2. Optical Spectroscopy 

In order to investigate correlations between behavior 
of CO lines, the photosphere, and the circumstellar com- 
ponents observed in atomic lines, we complement Table 
2] with seven high-resolution optical cchcllc spectra ob- 
tained during our continuous monitoring program of the 
star in 2000-04 (Jy)- 

The spectrum of 2000 July 19 was observed during the 
deep outburst minimum with the Utrecht Echelle Spec- 



trograph on the 4.2m William Herschel Telescope at the 
Canary Islands (Spain). The nominal wavelength resolu- 
tion is R ~50,000, with continuum signal-to-noise ratios 
exceeding 100. For more detailed discussions about cal- 
ibration procedures and a thorough analysi s of the out- 
burst spectra we refer to lLobel et all l|200l . 

The spectrum of 2004 October 23 was observed to fol- 
low up our CGS4 NIR observations of September 2004 
that revealed strong CO emission. It was obtained 5 
weeks later, and about a week before the CorMASS NIR 
observations, when the remarkable CO emission had al- 
ready started to decrease. The spectrum was obtained 
with the KPNO 4m Mayall telescope using the Echelle 
Spectrograph configured to provide R ~40,000 in the red 
(6000-8500 A in 22 orders). The wavelength calibration 
was performed with the Th-Ar lamp spectra, with a small 
additional shift in the zero-point as inferred from terres- 
trial O2 lines. 

The other optical echelle spectra of 2003-04 were ob- 
served with the Sofin spectrograph on the 2.6m Nordic 
Optical Telescope (La Palma, Spain). The spectra are 
observed with the medium-resolution Camera 2 (R ~ 
80,000) with a wavelength range covering about 5000 A 
over 42 echelle orders. Each observation combines a num- 
ber of subsequent exposures with two central wavelength 
settings until a minium S/N ratio of ~50 is reached. The 
S/N ratios for the echelle orders centered on the CCD 
can therefore exceed 300 in the wavelength region of the 
spectral lines of interest. Typical effective exposure times 
are limited to 10 min. per exposure to avoid CCD sat- 
uration for this optically bright star. Note for example 
that the large brightness of p Cas allows us to integrate 
high-quality optical spectra during dusk or dawn (i.e. at 
the beginning or end of a bright night run) with large 
zenith distances (>50 deg). For some of our spectra it 
yields large air-masses (e.g. ~1.6 for 2004 December 25) 
with strongly contaminating telluric lines, but which oth- 
erwise do not influence our spectral regions of interest. 
The Sofin spectra are calibrated with Th-Ar lamp spec- 
tra observed after each science exposure. The standard 
echelle reduction steps are performed using the 4A data 
acqui s ition and reduction software package described in 
lllvinl i|2000l) . Special consideration is given to the re- 
moval of scattered light inside the spectrograph, and to 
the extraction of spectral orders including the removal of 
cosmic spikes. The wavelength dispersion solutions yield 
mean accuracies better than ±1 kms -1 . The echelle 
spectra are not calibrated in an absolute flux scale, but 
special care is taken to apply the required blaze correc- 
tions for an accurate placement of continuum flux levels 
for continuum normalization. The local continuum lev- 
els are either determined with a polynomial fit through 
pre-determined anchor points in wavelength regions that 
are relatively free of spectral lines, or are based on a best 
visual estimate systematically applied to all of our obser- 
vations. Note that for these large S/N ratio optical spec- 
tra the accuracy of the placement of stellar continuum 
level is better than 1-2%, which is sufficient for measur- 
ing the small intensity changes we observe over time in 
weak absorption and emission lines f H5.1fl . 

3. REMARKABLE VARIATIONS OF p CAS NEAR-IR 
SPECTRUM 
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3.1. Spectra near Minimum and Maximum Light 

Due to a persistent lack of systematic libraries in 
the literature of NIR spectra of warm supergiant stars, 
we observed CorMASS spectra of several bright super- 
giants for comparison with p Cas. The spectra were 
obtained in 2003-04 with a setup identical to that for 
p Cas. Figure ^ shows our low-resolution NIR spec- 
tra of p Cas with less-luminous supergiants of similar 
spectral type: the non-variable p Per (GO lb) and the 
3d-period classical Cepheid Y Aur (F-G) in two vari- 
ability phases. HR 8752 is another semi-regular hyper- 
giant which was the spectroscopic twin of p Cas in the 
1970s but whose T s ff increased to abo ve 8000 K since 
then (Israclian, Lo bel fc Schmidt! 119991) . Due to the low 
resolution of CorMass it is difficult to determine if HR 
8752 currently is an early A or F star, although it clearly 
reveals a phase with strong emission observed in Ha, 
Pa/3 and the Ca II IR triplet lines. The latter confirms 
the strong emission we observe in the violet wings of the 
IR Ca II lines in 1998 October with NOT-Sofin, although 
these lines also reveal a central strong absorption com- 
ponent in the higher resolution spectra. The Ha line in 
HR 8752 has a permanent triple-peaked emission profile 
which, unlike in p Cas, may be attributed to the presence 
of a faint B-type companion discovered with IUE in HR 
8752. On the other hand however, the broad Ca II H & 
K resonance lines of both hypergiants lack strong cen- 
tral emission cores, signaling the absence of a classical 
chromosphere in these cool hypergiants. The more de- 
tailed comparison between the two hypergiants is given 
in fljTTl 

The overall spectral energy distribution of p Cas be- 
yond 1 pm compares well with what is expected for a 
F-G -type star (basically a black-body of T ~ 6000-7000 
K). It confirms the sm all reddening of only Ay=l m .3 
ijZsoldos fe Percvll|199Hn . of interstellar origin. The red- 
dening is barely noticeable with a slight flattening of the 
spectrum shortward of 0.8 pm. There is no substantial 
NIR excess around 2 pm as could be expected from warm 
circumstellar dust near the sublimation temperature of 
^1500 K. The low surface gravities (log g) of p Cas and 
HR 8752 are evidenced by strong O I 0.777 pm triplet 
lines in the low-resolution spectra. Other enhanced lines 
observed in p Cas around 1 pm belong to high-excitation 
transitions of C I and S I, and to low-excitation transi- 
tions of Sr II (the latter is also an atmospheric gravity 
indicator). These lines can signal high excitation condi- 
tions due to the low gas densities of hypergiant atmo- 
spheres, as well as a possible metal richness for p Cas (A 
log Z=0A ± 0.3; ILobel et alJ lll99l 'l. 

The K band spectrum of p Cas also shows CO absorp- 
tion. CO NIR bands usually occur only in stars of K- 
type or later l|Kleinmann fe Hall 1986: Wallace & Hinklc 
1996, 1997) and are indeed not detected in any other yel- 
low supergiant of our sample. Even more astonishing is 
the presence of strong CO band emission in 1998 Octo- 
be r (CGS4 archival s pectrum) when the optical spectra 
of ILobel et al.l l)2003f) reveal atomic emission lines. The 
CO emission can hence be considered as the continuation 
of p Cas' emission line spectrum from optical to infrared 
wavelengths, prompting a more detailed search for other 
NIR emission lines besides CO and Na I. 

In 2004 July we observe again a number of emission 



lines in the optical spectrum. We also observe the emis- 
sion lines, including CO, one month later with a medium- 
resolution SpeX observation. The most obvious emission 
features are marked in Fig. [3J We identify them as due 
to a low excitation energy transition (below 3 eV) of Ti I 
at 1.04 pm, while the lines at 1.5 and 1.7 pm of Mg I 
and of Na I at 2.2 pm emerge from considerably higher 
energy levels above 3 eV in the atom. NIR emission lines 
are more easily detected in F-G stars toward longer wave- 
lengths because the continuum flux is weaker and is de- 
void of strong absorption features. We indeed observe a 
number of emission lines in the L band. They are marked 
in Fig. |2Jand wavelengths and fluxes are listed in Tabled 
These lines cannot be identified with molecular features 
observed in cool giants, such as SiO, HC1, OH , or HCS 
l|Ridgwav et aUll984t IWallace fc Hinkkl l2002). Because 
of the sharp line profiles we think they are of atomic ori- 
gin. We searched for candidate transitions with atomic 
line properties comparable to the emission lines we al- 
ready identified in the optical and NIR spectrum. We 
queried the Atomic Line List database compiled by van 
Hoof (URL: www.pa.uky. edu/^peter/atomic) but ob- 
tain for each emission line more than one candidate 
atomic transition. One resolution element in our L-band 
spectrum equals about four times the average line width 
of emission lines (FWHM < 30 kms -1 as measured from 
the optical and NIR atomic emission lines). We there- 
fore refrain from further NIR emission line identifications 
until a higher resolution spectrum becomes available. 

We observe that the NIR emission spectrum appears at 
least twice simultaneously with optical emission phases 
that also coincide with phases of maximum light. Next 
we discuss the NIR spectrum for a phase of minimum 
light, which also occurs during the unusual outburst 
epochs. Figure [3| shows a NIR UKIRT spectrum of p 
Cas during the outburst in 2000 October. It is observed 
in the deep visual brightness minimum when the star was 
~1 magnitude fainter. To our knowledge it is the first 
IR spectrum of p Cas observed during a major outburst 
event. Compared to our SpeX spectrum observed around 
maximum brightness in September 2004, emission lines 
are not observed in outburst except for Na I 2.2 pm, 
whereas the absorption spectrum is much stronger (ex- 
cept for Br7) . The Na I emission lines protrude through 
underlying absorption lines. The spectrum is overall 
compatible with the photospheric spectrum of an early 
K-type supergiant. The central depths of CO, Ca I, and 
Al I lines, and of the Na I absorption lines (and possibly 
also a CN band but which is in the noisy region of telluric 
absorption) compare well with the K4-type supergiant 63 
Cyg in Fig. [21 This is in agreement with a T c g estimate 
of ~4250 K obtain ed for the optical outburst spectrum in 
ILobel et al.l l)2003fl . The outburst spectr um reveals how- 
ever a number of interesting peculiar properties. Br7 is 
unusually strong for a K-type star with a central line 
depth only slightly weaker than in the SpeX spectrum 
of T e ff ~7200 K. On the other hand, Ha is considerably 
weaker during outburst ijLobel et al.l 12003^1 . The low- 
excitation (xiow <2 eV) Ti I absorption lines at 2.17- 
2.23 pm are also very remarkable. They appear some- 
what stronger than in K- or G-type stars, in particular 
the 2.178 pm line that compares better to the M2-type 
supergiant p Cep with T D g ~3500 K. It is of note that 
these lines have been observed in emission in p Cas dur- 
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ing a shallow brightness minimum of 1980 June (Sheffer 
[Hal- They were also observed in emission for a short 
period of about one month during the eruption of the 
famous light-echo peculiar supergiant V838 Mon, while 
the object assume d minimum brightness and NIR CO 
was i n absorption ( Ban eriee & A shok 2002 iLvnch et alJ 
120041) . We therefore propose that the strong Ti I lines in 
the NIR UKIRT spectrum of p Cas can originate from 
a 3000-4000 K g as envelope ejecte d during the strong 
2000-01 outburst (|Lobel et al.ll2003|) . 

3.2. Two- component CO Lines 

We investigate the possible correlation of the NIR CO 
variability with the brightness curve utilizing all spec- 
tra available to us. The visual light curve in Fig. 0] 
is based on observations by the American Association 
of Variable Star Observers (AAVSO). Although their 
dataset is less precise than other photo-electric obser- 
vations, it is the only one available that continuously 
monitors p Cas from its earliest spectroscopic observa- 
tions until our most recent observations of 2004. We 
utilize the visual data to roughly estimate the variability 
phase (maximum vs. minimum light, brightness increase 
vs. decrease) as Srd variables are not strictly periodic, 
but only 'quasi-periodic'. Note that the light curve of 
this pulsatin g^ star follow s the radial movements of its 
photosphere. ILobel et al.1 (|2003f) showed that the visual 
brightness changes mimic the radial velocity curve of p 
Cas after 1997, with an offset of about a quarter of a 
period. The atmosphere begins to expand near maxi- 
mum light and reverses to contract after minimum light, 
for example as is also observed in more regular pulsating 
stars such as classical Cepheids. Unfortunately, radial 
velocity measurements require high-resolution spectro- 
scopic monitoring which is not available for all epochs of 
our observations. We therefore use the visual light curve 
as a proxy for determining the corresponding pulsation 
phases. 

Figure 0] plots the low-resolution K band spectra of 
p Cas against the light curve. We find that the CO 
band always appears in emission near phases of max- 
imum light (when the stellar radius already begins to 
expand) and that the band returns to absorption soon 
after the star dims. The CO absorption in the outburst 
spectrum of 2000 (shown in detail in Fig. [3Jl is only 
moderately strong, which is somewhat surprising consid- 
ering the deep brightness minimum. However, at the low 
spectral resolution one can not dismiss the possibility of 
filling-in by stronger than usual (residual) CO emission. 
Alternatively, the weakness of the CO absorption can be 
explained if the non-photospheric upper layers have been 
perturbed, if not entirely disrupted by the eruption. 

We further concentrate on the recent emission epoch 
of 2004 with the NIR echelle spectra. The emission 
epoch started around 2004 May (based on optical emis- 
sion spectra l|Lobel et al.ll2005(i ) and ended around Oc- 
tober. The resolution of the SpeX spectrum is too low 
to resolve individual line profiles of p Cas, although it 
provides a very large wavelength coverage. We therefore 
obtained a CGS4 echelle spectrum, with R^8 kms -1 per 
two pixels, covering five spectral regions centered on Bi'7, 
Na I 2.2 /im, and the first overtone Av 2-0 band of CO, 
from the bandhcad at 2.29 /im to 2.322 pm, where the 
band begins to overlap with the Av 3-1 CO band. The 



SpeX and CGS4 echelle observations are almost simulta- 
neous (ten days apart). 

The top panel of Fig. [S] shows our NIR echelle spec- 
trum (solid line) toget h er wi th a FTS spectrum from 
ILambert. Hinke fc Hall 1(1981]) (dotted line). The FTS 
spectrum was obtained during the brightness minimum 
of 1979 and reveals the CO band almost completely in 
absorption. The echelle spectrum, observed on the de- 
scending branch of the light curve before brightness min- 
imum, clearly harbors an absorption component in CO 
as well, although the band is dominated by emission in 
the low-resolution spectrum. 

The detailed profiles of several CO lines in 2004 are 
shown in Fig. [SJs, together with the Na I doublet and 
Br7 lines. The prominent Na I emission lines are very 
symmetric and centered around the center-of-mass ve- 
locity of the star. The Br7 line is in absorption and its 
core appears somewhat blue-shifted. The broader wings 
of this strong line extend to ^100 kms -1 in both direc- 
tions, also observed in Ha. The CO lines however consist 
of both emission and absorption components, with fluxes 
that systematically vary with the excitation energy of 
each transition. The low-excitation CO lines (with small 
Js) are in strong absorption, while high-excitation lines 
(with large Js around the bandhead) are in pure emis- 
sion. The CO emission line maxima appear to shift red- 
wards with decreasing Js. We think however that the CO 
emission lines are always centered around the star veloc- 
ity (or slightly blue-shifted by <7 kms -1 ) and the pro- 
gressive emission line redshift away from the bandhcad 
is only apparent due to a strengthening of the nearby 
CO absorption line components. With a FWHM of < 
25 kms , the CO emission component is narrower than 
the Na I emission lines and appears round-topped. The 
absorption component of CO is centered around Kadheiio 
= —80 kms -1 (or —33 kms -1 relative to the star veloc- 
ity), and has a prominent violet line wing extending up 
to at least —100 kms -1 (or —53 kms -1 relative to the 
star velocity). 

A detailed comparison of the two CO spectra observed 
25 years apart reveals a striking transformation from 
deep absorption lines in 1979 to apparent P Cygni pro- 
files in 2004. As is indicated by our low-resolution spec- 
tra, the transformation from absorption to emission ac- 
tually occurs on a regular basis every pulsation cycle. 
Sheffer's long-term high-resolution NIR monitoring of p 
Cas offers a unique opportunity to investigate its causes 
by tracking the detailed evolution of the CO line profiles. 
Figure [Hp shows a wide diversity in line shapes exhibited 
by CO in 1979-2004. For a fu l l set o f profiles we refer to 
Figs. 5.15 and 5.16 of lShefferl (jl£93). When we correlate 
them with the corresponding pulsation phases and com- 
pare with optical line profiles, a coherent interpretation 
can however be proposed. 

First, the variations of both emission and absorption 
line fluxes are larger when the amplitude of the light 
curve (A V) is also larger. Compare for example the 
strong CO lines observed in 1979-81 and 2000-04 with 
the 1987-89 epoch when the CO absorption lines remain 
permanently split. During the latter epoch pulsations 
are weak (see Fig. and the two absorption line com- 
ponents vary only slightly in intensity and Doppler shifts. 
Secondly, prominent CO line emission only occurs near 
phases of maximum light. It is always centered within 
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±10 krns" 1 around the star velocity. On the other hand, 
the CO line absorption can both be either red- or blue- 
shifted, with shifts of <30 kms -1 . The precise depen- 
dence of absorption on the corresponding pulsation phase 
however remains unclear. This is in part because the 
CO line absorption obviously becomes filled in with vari- 
able emission. The outer atmosphere of p Cas is also 
velocity stratified and only loosely coupled to the pulsa- 
tions of its lower photosphere, which can yield different 
Doppler displacements for the CO absorption and emis- 
sion line formation regions. Finally, it is important to 
point out that the properties and the time-scale of these 
CO line profile variations closely resemble the profile evo- 
lution observed in optical at omic split lines suc h as the 
Fe I 5506 A line in Fig. 2 of lLobel et all <|2003j) . When 
one of the split absorption components strengthens the 
other one weakens. The double CO core profiles can thus 
simply be interpreted as one broad absorption line with 
semi-regular Doppler displacements, in which core split- 
ting occurs when it passes across a static central emission 
reversal inside the line core. 

The broad range of NIR CO line morphologies in p Cas 
argues firmly against a P Cygni line profile formation re- 
gion. The CO emission flux varies on time scales directly 
linked to the variability cycles of the photospheric aborp- 
tion spectrum. This is unlike classical P Cygni profiles 
observed in ultraviolet resonance lines of hot stars that 
are static and almost invariable in time. It signals that 
the major portion of the CO emission line formation re- 
gion is influenced by the recurrent photospheric changes 
in p Cas, rather than for example emission emerging from 
a very far extended optically thin circumstellar enve- 
lope that produces redshifted photons by backscattering 
in the far-side expanding hemisphere. The 2.3 pm CO 
emission can therefore not emerge from a (possibly) very 
extended molecular (e.g. CO) gas envelope around the 
supergiant. This is further evidenced by NIR CO line 
profiles o bserved i n HR 8752, another variable cool hy- 
pergiant. iShefferl (|1993|) observed that either only emis- 
sion occurs in its CO lines, or that it is accompanied by 
red-shifted abso rption, yielding an inverse P Cj/qm-type 
profile instead <|Lambert. Hinke fc HaiH ll981I The 2.3 
pm CO feature however vanished in HR 8752 in the late 
1980s, approximately at the same time when its optical 
variability became much smaller and the effective tem- 
perature increased f< J5.2.1[l . We now investigate where 
the CO emission line formation region is located in p 
Cas. 

4. MODELING CO EMISSION 

Since the 2.3 pin CO band is not observed in 'nor- 
mal' F-G -type supergiants, neither the CO emission nor 
the absorption components appear to be of photospheric 
origin in p Cas. The supergiants of our comparison sam- 
ple however have smaller luminosities with larger gravity, 
and molecular bands can be strongly gravity-dep endent 
ijKleinmann fc Hal]|ll986l lMcGovern'etaIll20o|) . Fur- 
thermore, since p Cas is a pulsating star, can the 
strong variability of the underlying photospheric spec- 
trum cause apparent variability in CO, as well as in other 
'external' spectral lines? The strong correlation of the 
CO emission flux with variability phase signals its for- 
mation close to the pulsating photosphere. We therefore 
adapt a numerical model of the average photosphere to 
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probe a number of spectral properties of the CO band 
observed in the high-resolution CGS4 data. 

4.1. Photospheric Parameters 

Every pulsation phase of p Cas requires an estimate 
of at least two photospheric model parameters: T e g and 
log g. We first compare the CorMASS spectra of p Cas 
with its less-luminous cousins, the Cepheids and the non- 
variable yellow supergiants of our sample. In the latter 
stars the photospheric lines are sufficiently narrow to ap- 
ply the classical method of absorption line depth ratios to 
reliably infer the T c g. Figure ^ shows two spectra of the 
3d-period classical Cepheid Y Aur obtained about half 
a pulsation cycle apart. High-resolution optical spectro- 
scopic studies show that T e g of short-period Cepheids 
can yary between 5500 K and 6 500 K with log g= 1.9- 
2.3 l|Kovtvukh fc Gorloval feOOOj) . For the early G non- 
variable p Per in Fig. Q]we adopt T c g = 5330 K and log 
g = 1.8 (V. Kovtyukh, private comm.) When comparing 
CorMASS spectra of p Cas of 2003 November and 2004 
October to Y Aur and p Per we find that T c g should ex- 
ceed 5500 K for p Cas. To get a better constraint on the 
temperature, we fit the CorMA SS spectra of p Cas (a f- 
ter dereddening by A v = l m .3 (|Zsoldos fc PercvllT99Hn 
with the Next-Ge n -giant model spectra c omputed with 
the Phoenix code l|Hauschildt et al] ll999h Our best fits 
confirm a log g value of ^0.0, whereas for T e g only a 
lower limit can be derived due to a lack of models with 
Teff > 6600 K and log g < 2.0 (Fig. EJ. We obtain a 
T eff > 6400 K for the spectrum of 2003 November, but 
slightly cooler for 2004 October. The light curve (Fig. 
0J) shows that the star assumes maximum brightness be- 
tween 2003 November and 2004 October. The star is 
dimmer in the V band in 2004 October compared to 
2003 November, which confirms the smaller T c g derived 
from our best fit procedure. If one assumes that the CO 
absorption component forms in the photosphere and be- 
comes stronger with smaller T e g, then the shallower CO 
absorption band observed in 2004 October must result 
from filling in by emission which occurs on the descend- 
ing branch of the light curve. In the following Section we 
investigate whether or not the CO absorption component 
forms in the stellar photosphere. 

We can als o constrain T R g from our high-resolution op- 
tical spectra. lLobel et alJ l)2003j) established a tight rela- 
tionship between the equivalent width (EW) of the Fe I 
5572 A absorption line and the T e g of p Cas. Our optical 
echelle spectrum of 2004 April 6 indicates an EW for the 
Fe I line of 348 mA which corresponds to T off =7290 K. In 
the spectrum of 2004 October 23 we measure an EW=465 
mA yielding a T e g =6910 K. In addition, the visual 
brightness in 2004 October was comparable to 2003 June, 
for which Fe I provides ^6800 K. These estimates agree 
well with the lower limits estimated from the CorMASS 
spectra. They confirm that the hypergiant was cooling 
down while dimming to minimum brightness (after early 
spring 2004), but it did not cool to below 6000 K. We can 
hence adopt a reliable estimate of T c g =7200±100 K for 
the SpeX and CGS4 emission spectra of 2004 Septem- 
ber, and of 6800±100 K for the CorMASS absorption 
spectrum of 2004 October. 

4.2. Temperature Minimum in the Upper Atmosphere? 
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We use a Kurucz model of the photosphere with 6750 
K and log 3=0.5 to compute the first overtone absorp- 
tion band of CO. We adopt a value of 0.5 because it is 
the smallest gravity model available for this temperature. 
We perform LTE radiative transfer syntheses of the 2.29- 
2.32 /im wavelength region using Kurucz line lists. The 
input line data include the 2.3 pm CO band transitions 
and a number of strong atomic transitions in this NIR 
region. The top panel of Fig. [7|shows the kinetic gas tem- 
perature structure of the model. The bottom panel plots 
the computed spectrum (dot-long- dashed line). The com- 
puted spectrum only shows some weak atomic absorption 
lines, without discernible CO features. We checked that 
noticeable CO absorption can only be computed with 
models of T e g < 6000 K, significantly below the values 
we determine from our spectral monitoring in 2003-04. 
These hydrostatic models of the photosphere therefore 
do not produce strong CO absorption bands in F-type 
stars such as p Cas. The kinetic gas temperatures in the 
photosphere are too large to yield a sufficient CO abun- 
dance, even in the low-gravity conditions of super- and 
hypergiant atmospheres. 

Since we observe that the CO emission variability cor- 
relates more strongly with the stellar pulsations than 
the CO absorption, we model the emission with atmo- 
spheric layers in the vicinity of the photosphere. We 
adapt a Kurucz model for the 2004 September emission 
spectrum with T off =7250 K and log g=0.5 (Fig. |7| dot- 
short-dashed line) by adding a kinetic gas temperature 
minimum in the upper photospheric layers over a small 
range of optical depths (long-, short-dashed and dotted 
lines). The modified model produces CO emission when 
the Tg-minimum is placed below log m < —4 in the col- 
umn mass scale. The observed emission flux at the CO 
bandhead does not appear to be influenced by absorp- 
tion and can correctly be computed if the Tg-minimum 
is located around layers with log m = — 3. 1±0. 1 . The lat- 
ter column mass value corresponds to only ~5% of the 
mean stellar radius above the photosphere. We apply 
a projected microturbulence velocity broadening of 11 
kms -1 to compute the NIR spectrum. The supersonic 
microturbulence velocity fits the observed CO emission 
line widths, without invoking extra macroturbulence or 
rotational line broadening. Our composite model pro- 
duces the observed CO emission flux. It however fails to 
precisely match the detailed shape of the complete CO 
band. For example, the computed spectrum obviously 
lacks sufficient flux longward of the bandhead. Further 
improvements of the fit would require an even more de- 
tailed atmospheric model with additional modifications 
in the other thermodynamic conditions besides the T s - 
structure (in the current model they are only roughly 
interpolated to match the extension of the upper temper- 
ature profile). More sophisticated and detailed non-LTE 
radiative transfer calculations of CO would be appropri- 
ate in the very low density conditions of our model, but 
they are beyond the scope of the present investigation. 

How can the temperature minimum required in these 
model fits be produced in the outer atmosphere of p Cas? 
A similar atmospheric temperature structure has been 
proposed for the Sun to explain the fundamenta l band 
of CO at 4.7 pm observed off-limb in emission l)Avresl 
2002). The solar CO emission forms in the so-called 
'CO-mosphere', a transitional region located between the 



photosphere and the chromosphere at an altitude of 600- 
1000 km. This intermediate emission region is an inho- 
mogeneous gas mixture with an average T g of ~3000 K 
permeated by chromospheric 'fingers' of T g ~10,000 K. 
There is however no evidence of a classical chromosphere 
in p Cas: the most important stellar chromospheric indi- 
cators such as emission in the cores of the Ca n H & K 
and Mg n h & k lines have never been observed. 

Because we observe the strong NIR CO emission only 
during variability phases of mean atmospheric expan- 
sion it may indicate an alternative dynamical mecha- 
nism that can cause temporal sharp structures in the 
temperature- and density-profile of the upper atmo- 
sphere. An outwardly propagating field of weak stochas- 
tic shock wave trains was for e xample proposed by 
Ide .lager. Lobel & Israelianl l|1997fl to model the super- 
sonic microturbulence velocities observed in the extended 
atmosphere of p Cas. The CO emission could be pro- 
duced by a pulsation-driven circumstellar shock wave 
in the cooling layers behind the shock front where the 
shocked gas expands while flowing into the following 
shock wave. After shock passage the cool CO layers 
are heated by compression in the trailing shock which 
partially dissociates the molecular CO gas fraction and 
quenches its local emission. 

5. DISCUSSION 
5.1. Comparing CO with Peculiar Atomic Lines 
5.1.1. Atomic Line Classes 

To compare the NIR CO lines with atomic line profiles, 
we plot in Fig. [8] representative examples for the main 
line morphology groups observed in p Cas. The NIR 
lines are from the UKIRT spectrum of 2004 September 
17 with emission lines during a phase of declining bright- 
ness. The optical lines are from the KPNO spectrum 
of 2004 October 23 when the emission lines weaken be- 
fore brightness minimum. Both spectra have comparable 
spectral resolution of ~8 kms -1 , sufficient to resolve in- 
dividual line profiles. Although the emission spectrum 
weakens in 2004 October (see the CorMASS spectrum in 
Fig. and the atmospheric expansion must have decel- 
erated somewhat, we believe that the overall line shapes 
and Doppler shifts do not appreciably change between 
the two observations (less than ~1/12 of P apart). We 
can therefore directly compare the NIR and optical line 
shapes, within the same pulsation cycle and nearly the 
same phase, which is crucial for this long-period semi- 
regular variable star. 

In Fig. [S]we overplot spectral features of p Cas with 
high-resolution spectra of Arcturus (a Boo) and the Sun 
(both aligned to the rest velocity of p Cas), because their 
spectral lines are narrower and more symmetrical. These 
reference spectra permit accurate line identifications in 
p Cas whose lines are very broad and often blended or 
distorted with emission. Furthermore, the T e g of ~4400 
K (a Boo) and 5800 K (Sun) provide a useful estimate of 
the sensitivity of the photospheric absorption spectrum 
with the T e ff changes of p Cas, which can reveal the in- 
trinsic variability of emission lines in different pulsation 
phases. We compare with Arcturus and the Sun despite 
the larger atmospheric gravity values of log g— 1.5-2.0 
ijBell. Edvardsson. fc Gustafssonlll985|) and 4.4, respec- 
tively, than in p Cas (log g ~ 0) because detailed syn- 
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thetic spectra of cool hypergiants are presently not avail- 
able. Although there is no clear relation between all the 
spectral lines of these three cool stars (e.g. strong absorp- 
tion lines in the Sun and Arcturus are sometimes absent 
in p Cas), the comparison reveals that the diverse line 
profile shapes in p Cas are intrinsic to the peculiar spec- 
trum of this hypergiant, and do not result from common 
absorption line blending. 

The optical lines of p Cas can be divided into four main 
types we discuss in more detail: 

I. Photospheric absorption lines of neutral and singly 
ionized metals (e. g. Fe I 5572 A used as T ff indica- 
tor in H4.1|l . Hydrogen absorption lines can be included 
in this group as well, although their line formation re- 
gions extend far above the photosphere, while Ha can 
ex hibit emission co mponents in both line wings (Fig. 2 
of lLobel et afll200l . 

II. 'Split' absorption lines with a central emission peak 
that is essentially constant in both flux and wavelength. 
In the stronger lines of multiplet (with larger oscillator 
strengths log(gf)) the wings of the emission core show 
up against the broad and intensely saturated underly- 
ing absorption cores (example, Fe I 5269 A line). The 
cores of weaker lines (such as Fe I 6191 A) split only 
when the absorption core crosses over the central emis- 
sion peak, which occurs when the radial velocity of the 
photosphere is similar to its mean. The FWHM of unsat- 
urated split lines is comparable to single absorption lines 
(^60 kms -1 ) and indicates a common photospheric ori- 
gin for both. 

III. Metallic emission lines of variable strength. They 
appear prominently above the level of the stellar contin- 
uum soon after maximum brightness, but they diminish 
and then completely vanish around minimum light. A 
good example of this type of emission line is the semi- 
forbidden Fe I triplet around 0.8 /an. 

IV. Metallic emission lines of roughly constant 
strength. This class includes a few [Ca n] lines and the 
2.2 fim Na I emission lines. 

Figure |H1 shows a number of illustrative lines for each 
category. They represent a range of line strengths and 
are minimally blended. The lines progressively change 
strength from absorption to emission rather than repre- 
senting separate line groups. The detailed shape of each 
line appears to depend on the excitation energy, and to a 
lesser degree on the oscillator strength. Disregarding all 
hydrogen lines and considering only lines with xiow below 
a few eV, we find that the emission component strength- 
ens while the absorption component weakens with de- 
creasing log(gf). The Na I 2.2 pm emission lines are 
the only exceptions because these moderate excitation 
lines do not show photos pheric absorption in F-G stars 
ijKleinmann fc Halilll986l our Fig. UJ. 

5.1.2. Putting CO Lines in the Context 

The NIR CO lines contain both emission and absorp- 
tion components with very similar profiles compared to 
the type II and III atomic lines. Note however that it 
has not been established yet if there is any physical dif- 
ference between the type II and III lines, if they form 
in two separate atmospheric regions or even with differ- 
ent excitation mechanisms. The morphologic similarity 
of both line types is apparent with a smooth transfor- 
mation from emission below the local continuum level in 



the strongest split lines to (prominent) emission above 
the continuum for the weaker type III lines (Fig. |HJ. 
For example, the central emission cores in the split Na 
D lines can only rise a bove the continuu m level during 
the outbursts of p Cas llLobel et al.l l2003). Other exam- 
ples of these mixe d line types are presented in Fi g. 2 of 
lLobel et al.l l)2005|) . Based on the line similarities. lLobel 
l)1997|) proposed that emission components in both types 
emerge from two circumstellar conic shells formed at the 
interface of a supersonically expanding bipolar wind that 
collides with previously ejected material from the hyper- 
giant's violent mass-loss history. At the interface the 
wind velocity decelerates, becomes subsonic, and forms a 
steady (or standing) shock wave in the outer atmosphere. 
A curve of growth analysis applied to optical multiplets 
of the central emission cores in split Fe I lines yields an 
excitation temperature of 3050 K for the shocked region. 
Provided that these central emission lines are sufficiently 
optically thin, classical shock wave theory yields a thick- 
ness for the bipolar conic shells smaller or of the order of 
the mean stellar radius (using an average mass-loss rate 
of M=9TO _5 M0 derived from the far violet extended 
wings of photospheric absorption lines). The observed 
emission fluxes yield distances of 15 to 150 i?* f rom p 
Cas f or the location of the emitting shock interface l)Lobel 
119971 p. 54) for realistic cone opening angles of more than 
30 deg. Smaller opening angles for bipolar outflow would 
provide distances to the shock interface that exceed es- 
timates in the literature fo r the distance to the dust for- 
mation region near p Cas l).Iura & Kleinman"nlll990ri . 

At these large distances, strong influences of the pul- 
sating photosphere on the emission excitation region are 
not expected, consistent with the constant flux of the 
central emission cores of split lines (type II). The Doppler 
stability of the emission lines with pulsation phase results 
from the subsonic outflow velocity behind the standing 
shock wave. We however observe different properties for 
the type III emission lines. These emission lines can ap- 
pear sporadically in the core of unsaturated absorption 
lines and their shapes can alter while the emission varies 
with the pulsations. It is important to know if these 
lines are intrinsically variable in total flux and Doppler 
position. The top panel of Fig. shows the time evo- 
lution of the unblended emission-free photospheric Fe I 
5572 A line during a complete pulsation cycle from 2003 
June shortly before minimum light, throughout the max- 
imum of 2004 May, until the minimum of 2005 January. 
The outburst spectrum of mid 2000 is also shown. The 
violet line wing develops at the beginning of the bright- 
ness decline (the atmospheric expansion can be traced 
by the blue-shifted line centroid) and shrinks back when 
the star reaches minimum light. The EW (and central 
depth) of the line increases monotonically from maxi- 
mum to minimum light because T e ff decreases. The bot- 
tom panel of Fig. ED shows a spectral region with type 
III lines in the same pulsation cycle. While the absorp- 
tion spectrum steadily strengthens from 2004 April to 
December (e.g. the Fe I 5572 A and the long-wavelength 
wing of Eu II 6645 A become stronger), the emission 
component clearly brightens as well between April and 
July and fades by December. The emission line evolu- 
tion resembles the blue wing behavior in the top panel 
of Fig. [5J The variations of the absorption line wing 
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trace changes of the wind opacity which also enhances 
between April and October 2004. The prominent type 
III emission lines that peak above the continuum there- 
fore vary intrinsically. Their variability is linked to the 
photospheric pulsations, since the emission flux increases 
together with the wind opacity during phases of fast at- 
mospheric expansion. 

The detailed evolution of the NIR CO resembles more 
the type III than the type II lines. The CO emission flux 
tightly correlates with the pulsations, signaling it is di- 
rectly influenced by the photospheric dynamics indicated 
with our models in 21 On the other hand, type II emis- 
sion lines appear invariable and must form further away 
from the surface. There are also slight morphological dif- 
ferences between type III/CO and type II emission lines. 
The former lines are centered around stellar rest velocity, 
while the latter are blue-shifted by about —5 kms -1 due 
to the subsonic outflow behind the steady shock inter- 
face. The type III and the CO emission lines also appear 
more round-topped compared to the triangular shapes of 
the strong central emission cores in the optical split lines 
(Fig. |SJ. The former lines signal a spherically symmetric 
line formation region, while the latter lines can emerge 
from a more cone-like outflow observed nearly pole-on. 
The differences between the two line groups in Doppler 
shift and line shape are rather subtle, but may be crucial 
to determine whether or not they belong to two separate 
emission regions around p Cas. 

The variability of type III and NIR CO emission lines 
signals formation regions close to the pulsating photo- 
sphere of p Cas. But how close? Our model of the CO 
emission places the formation region at only ~5% of i?* 
because the gas density of the emission region decreases 
rapidly farther above the surface. Figure reveals how- 
ever that the optical emission flux increases together with 
the wind opacity determined from the far violet wings 
of the photospheric absorption lines. The properties of 
the stellar wind dynamics are determined by the photo- 
spheric pulsations and appear to be linked to the emis- 
sion line formation region as well. If one could observe a 
time delay between the increase of the violet absorption 
line wings (the mean wind opacity) and the onset of the 
emission line spectrum, it would directly provide the dis- 
tance the wind can travel between the stellar surface (e.g. 
the base of the stellar wind in the photosphere at log(r 
~2/3)) and the actual emission line region. A geomet- 
ric extension of the wind to ~2.5 i?* QLobel et alJH998j) 
has been modeled from the photospheric absorption lines 
of p Cas, which would disagree with the much smaller 
distance we determine for the NIR CO. Though, part 
of this difference may arise because our models for the 
CO band emission assume a LTE approximation which 
breaks down in the very extended atmosphere and wind 
of the hypergiant. If on the other hand one observes that 
emission spectrum develops earlier than the violet wings, 
then it would favor our shock hypothesis. To discrimi- 
nate between the two scenarios, a very precise timing 
of the detailed line profile evolution using more frequent 
spectroscopic monitoring intervals will be required. 

5.1.3. Na I IR Emission 

We close this section with a discussion of the Na I 
2.2 pm doublet lines. These NIR lines are always 
observed in emission in p Cas, even during the deep 



outburst minimum of 2000-01 (Fig. |3J). The two 
emission lines of the Na I 2.2 p,m doublet are of nearly 
equal strength. On the other hand, the spectra of 
non-emission supergiants in Fig. [3] show that the 
short-wavelength absorption line is stronger (due to 
larger log(gf) value and blending with other lines) than 
the longer-wavelength one, causing the short-wavelength 
emission component in the p Cas UKIRT outburst 
spectrum of 2000 October to appear less conspicuous 
at the medium resolution. The NIR Na I emission is 
also observed in many other cool la and IaO super- 
giant s, incl uding HR 8752 l |Laml^gr^. Hinke ^^igXt 
198ll iMcCregor. Hvland fc HillieH 119884 

Hrivnak. Kwok fc Geballd 11994 our Fig. [TJ. It is 
traditionally explained in the literature to result from 
pumping by stellar UV radiation, because of the high 
excitation energy compared to other optical emission 
lines. The excitation energy of the upper level is 3.8 
eV compared to ~2 eV for optical emission lines. The 
upper energy level for the Na D core emission can be 
populated by cascading from the lower level of the 
NIR Na I doublet through the 1.138 and 1.140 ^m 
transitions. Unfortunately, the latter lines fall in a 
wavelength region with strong telluric absorption and 
cannot reliably be identified in our SpeX spectrum. 
Calculations of fluorescence were originally carried out 
for a nother star - IRC+10420 l)Thomnson fc Borosonl 
Il977|) . This cool h ypergi a nt has a very large mass loss 
rate in quiescence (HEE3J- Thompson Borosonl lfl97l 
considered Mg I lines with relatively high excitation 
energy (xup = 6 eV) at 1.5 and 1.7 pm, and concluded 
based on their strength that they must be optically thick 
and therefore emerge from a region not far from the 
photosphere with sufficiently large electron density, such 
as a stellar chromosphere. There are however no clear 
indications of a classical chromosphere in p Cas ( H4.2|l . 
The Mg I lines in the H band of p Cas also appear much 
weaker compared to Na I, unlike in IRC+10420 (e.g. 
compare the SpeX spectrum in Fig. J H with Fig. 7 of 
iHumohrevs. Davidson fc Smith! l)2002f )V It is therefore 
not clear if one can postulate the existence of a classical 
chromosphere in all cool supergiants with the 2.2 pm. 
doublet observed in emission. The rather triangular 
shapes of the Na I 2.2 pirn emission lines can also point 
to a line formation region in a collimated outflow, 
rather than in a very extended spherical gas envelope 
around p Cas where the more round-topped forbidden 
[Ca n] lines form. Independent of the mechanism 
producing this permanent and strong line emission, the 
Na I 2.2 pm doublet serves as a powerful indicator for 
identifying extreme supergiants in the NIR, even from 
low-resolution spectra. 

5.2. Comparing CO Emission and Split Lines in p Cas 
with Other Cool Luminous Stars 

The CO absorption lines we observe in p Cas can orig- 
inate in the circumstellar gas envelope, whose existence 
has been inferred for cooler K-M supergiants from spec- 
tral modeling and inte r ferometric NIR m easurements 
(e.g. , iTsul Il988l l200lt IPerrin et all l200l ) . Near-IR 
molecular emission bands observed in p Cas are however 
a much rarer spectral phenomenon. The emission excita- 
tion mechanism is not likely caused by photo-ionization 
from a hot companion star because very long-term ra- 
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dial velocity monitoring does not reveal any indication of 
binarity, while recent far-UV spectroscopy with FUSE 
confirms that p Ca s has no significant UV excess either 
l|Lobel et al J 12005). as is for example observed for the 
hot companion in HR 8752 with IUE. Variable emission 
lines in cool stars however are often associated with wind 
variability, bulk movements of circumstellar gas due to 
wind outflow, or to periodically shocked layers observed 
in pulsating stars. CO Av 2- emission lines are ob- 
served in low-mass protos tars llBisca va et al .lll997l). in 
hot and cool supergiants ([Mc Gregor. Hvland . fc Hillicr 
11988atl:lKraus et alJj200ffl : iMozurkewich et all l|1987tk 
iMatsuura et alJ l|2002tl ). and in s ome extreme cases o f 
He-flash ob jects (like V838 Mon. iRushton et alll2005|) 
and novae {Ferland et alJll97ft iRudv et al.ll200^) . The 
wind and shock in p Cas however must be much more op- 
tically thin because of the observed absence of emission 
in the hydrogen lines and of prominent P Cygni-type line 
profiles (except for Ha). 

5.2.1. Cool Hypergiants 

How does the NIR spectrum of p Cas compare to other 
luminous cool hypergiants in the upper HR-diagram? 
The only two objects that at one time have spectro- 
scopically been comparable to p Cas are the hypergiants 
IRC+10420 and HR 8752. IRC+10420 is only ~0.1 dex 
more luminous in log(L), while HR 8 752 is less lumi - 
nous by the same amount than p Cas (jde Jaeerl ll998). 
Despite the close similarity of the basic stellar parame- 
ters, the three hypergiants differ in many other respects. 
IRC+10420 has a constant mass loss rate which is an or- 
der of magnitude larger than the other two, and is also 
a maser source of the earliest spectral type known. The 
NIR s pectra of 1976 llThomnson fc Bor oson 1977) and 
1984 l<Fix fc Cobb|ll987ft revealed Na I emission lines, 
with the higher members of Brackett series in absorp- 
tion, while Br7 was not detected. The weakness of 
the hydrogen lines at NIR wavelengths was explained 
with veiling by the spatially resolved dust envelope. In 
1992 lOudmaiier et al.l l)1994D discovered how blue-shifted 
emission features protruded through the NIR hydro- 
gen absorption lines, completely filling them in by 2000 
({Humphreys. Davidson fe Smith! r2002j) . Interestingly, a 
weak CO 2.3 pm emission band can be present in the 
latter spectrum, although it is observed at noise levels. 
The optical spectra show signs of both out- and down- 
flows. The var iability of this co o l hyp ergiant remains 
poorly studied. iCottlieb fc Lillerl l)1978fl present photo- 
graphic photometry signaling a strongly variable object 
before 1920 which was steadily increasing in brightness 
until at least the late 1970s. To reconcile all th ese obser- 
vations Humphreys. Davidson & Smith (2002) suggested 
that the variable spectrum of IRC+10420 could entirely 
originate from an opaque circumstellar wind rather than 
the photosphere, mimicking rapid evolutionary changes 
in the HR-diagram near the cool boundary of the Yellow 
Hypergiant Void. 

HR 8752 has a less extreme mass loss rate than 
IRC+10420, making it more interesting to compare 
with p Cas. It has its own peculiarities though. The 
optical and NIR spectra are long-term variable as in 
IRC+10420. The spectral changes however have dif- 
ferent properties. The emission spectrum has become 
systematically weaker while the absorption spectrum 



and colors indicate a steady increase in T e ff of ^1000 
K per decade over the past 30 years. The pho- 
tosphcric pulsations have almost ceased while the V 
amplitude decreased from ±0 m .2 in the early 1980s, 
to ±0 m .05 in the mid 1990s iPercv fc Zsoldol Il992t 
iNieuwenhuiizen fc de ,Iagerll2000t) . From 1960 until 1980 
it was in all aspects very similar to p Cas, includ- 
ing periodic line doubling and the occurrence of emis- 
sion lines (interestingly also in NIR CO), with both 
emission and absorp tion components that altered with 
the pulsation phase llHar mer, Lawson fc Sticklandl 1978 ; 
i Lambert fc Lucid fl97St ILambert. Hinke fc Halll Il981 : 
ISheffer fc Lambertill987t lShefferll!993F ~ 

We compar e the NIR CO spectra of HR 8752 from 
Sheffer (1993) with the brightness curve and observe that 
the CO emission lines appear at the stellar rest velocity 
while the emission flux peaks around maximum bright- 
ness, as we also observe in p Cas. The CO absorption 
varies less regularly and may occur both red- and blue- 
shifted. Surprisingly, by 1988 the NIR CO band disap- 
pears completely, followed by type II and a flux decline 
in the type III emission lines. The absorption lines re- 
main however very broad as observed in the optical spec- 
trum of HR 8752 in 1998 bv llsraelian. Lobel fc Schmidt! 
(1999). We observe no trace of NIR CO in our Cor- 
MASS spectra of 2003 November and 2004 October (Fig. 
nj. The Na I 2. 2 urn emission lines are still observed, 
although iShefferl l)1993D reported a decrease of the line 
fluxes. Unlike p Cas, HR 8752 has a main-sequence early 
B-type companion that may contribute to the hyper- 
giant's emission spectrum. For example, the [N n] 6548, 
6583 A lines are thought to originate in a Stromgren 
sphere surrounding the companion, but we do not ob- 
serve them in p Cas. What caused the disappearance of 
CO and other low-excitation lines observed before 1987 
in HR 8752? We suggest that either the decrease of the 
photospheric pulsation amplitude weakened the shock 
wave(s) that can excite these emission lines, or that CO 
molecules close to the photosphere dissociated with the 
large increase of T c g during following decades. But what 
can cause this very fast increase of T e s? Possibly we 
observe this hypergiant in transit across the Yellow Evo- 
lutionary Void, on its way of becoming a blue supergiant. 
Or perhaps its atmosphere is warming up with the ap- 
proach of the hot companion. Another possible expla- 
nation is that HR 8752 exhibits long-term secular vari- 
ability, perhaps comparable to the changes from more to 
less regular atmospheric oscillations we also observe for 
p Cas. Continued spectroscopic and photometric moni- 
toring of these unique luminous cool objects is needed to 
address these questions in further detail. 

5.2.2. Post-AGB Supergiants 

Hypergiants are very exceptional luminous stars. Fur- 
ther insights of their NIR spectroscopic variability can 
be gleaned from more abundant and more regularly 
pulsating less-luminous cool supergiants. Good exam- 
ples of stars for which we also observe NIR CO in 
emission are a number of yellow post-AGB stars and 
some proto-planetary nebulae precursors. It is un- 
clear at present if NIR CO emission occurs more fre- 
quently i n variable stars or for a certain spe ctral type or 
subclass. iHrivnak. Kwok fc Geballel l)1994|) found three 
proto-planetary nebular supergiants with CO emission. 
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The e arliest is GO la IRAS 22223+4327, a pu lsating vari- 
able. IMcOregor. Hvland. fc HilTiel (|1988bl) found CO 
emission in a few early-type Ia-supergiants in the Mag- 
ellanic Clouds, but they are not known to be photomet- 
rically variable. The earliest reported case is the RV 
Tau F4-type Ibp star AC Her (HD170756), which also 
shows Br7 in absorption and Na I 2.2 pm in emission 
similar to p Cas. NIR CO in this star switched from 
emission to absorption in three low-resolution spectra of 
lOudmaiier et al.l 0995). We have reduced an archival 
CGS4 spectrum of another post-AGB Srd variable star 
with NIR CO emission, HD 179821 (G5 la). The spec- 
trum was obtained on the same night of 2000 October 18 
as the outburst spectrum of p Cas. It shows the Na I 2.2 
pm. lines in emission but not a trace o f CO, confirming, 
the CO variability previously noticed bv lOudmaiier et aLl 
IIT9951 . 

5.2.3. R CrB Stars and p Cas Outbursts 

The NIR CO bands observed in R CrB stars offer an in- 
teresting comparison to p Cas. These less-luminous erup- 
tive variables have carbon-rich and hydrogen-deficient at- 
mospheres, but exhibit many similarities to p Cas other- 
wise. They cover the same range of T e ff, reveal promi- 
nent emission lines at times, and are known pulsators 
that m ost important ly exhibit deep visual brightness de- 
clines IIClavtonlll996() . It is thought that occultations by 
dense dusty circumstellar clouds cause the abrupt dim- 
ming. The dust clouds should occur rather close to the 
photosphere, dictated by the short time scales of the sud- 
den brightness decreases. But how can dust form at pho- 
tospheric g as temperatures above 5000 K in these F-G 
type stars? iHinkle. Jovce. fe Smithl l)1995j) observed the 
CO 2.3 ^m band in FG Sge, a post-AGB and He-flash 
supergiant, two months before it went into a deep bright- 
ness decline, which turned it into a typical R CrB star 
thereafter. They observe split CO line profiles with one 
absorption component near the stellar rest velocity and 
the other blue-shifted by 21 kms -1 . The CO lines are 
interpreted as evidence of a long-sought reservoir of cool 
material forming near the stellar surface. After the first 
circumstellar gas shell forms at a distance of ~2 R*, the 
second shell cools below 1000 K and forms dust while be- 
coming accelerated outwards by ra diation pressure on the 
dust. Based on these observations lGonzalez et all l|1998fl 
proposed that CO gas might be the leading coolant ini- 
tiating the nucleation of dust grains very close to the 
photosphere. H owever, a very rec e nt ext ended study of 
R CrB stars bv lTenenbaum et alJ l)2005f) did not estab- 
lish any clear correlation between the appearance of the 
2.3 pm CO band and phases of large brightness decrease. 
In fact, they model the observed CO band strengths with 
pure photospheric models. Their work points to a strong 
connection with p Cas for which we observe that the 
absorption line splitting is due to a pulsation phase de- 
pendent superposition of absorption and emission line 
components that recur on a regular basis, rather than be- 
ing direct tracers of cool circumstellar material expelled 
during outbursts. 

There are convincing indications that the trigger for 
dust formation i s actually locat e d in the photosphere 
of R CrB stars. IGonzalez etaD i|1998fl found that the 
photosphere cools by at least ~1000 K during the fast 
brightness declines. This atmospheric phenomenon is 



often neglected because of difficulties in interpreting R 
CrB spectra observed at brightness minima. The ab- 
sorption lines appear to weaken and to become dis- 
torted, possibly because the photospheric spectrum is 
blurred by the scattering from the expanding dust clouds 
(although these effects are not clearly observed in p 
Cas). The atmospheric cooling itself appears to result 
from some kind of dynamic t ransformation or instability. 
IKameswara Rao et alJ l)1999fl (see their Fig. 1) observe 
that at the beginning of a strong brightness decrease 
in R CrB the low-excitation photospheric lines reveal 
an emission component and become strongly red-shifted. 
The high-excitation lines in contrast become blue-shifted. 
This signals that the upper atmosphere decouples from 
the regular pulsation movements. The photosphere col- 
lapses, causing a stronger than usual compression, which 
in turn results in a stronger than usual atmospheric ex- 
pansion with a large brightness decrease. A very similar 
behavior is observed during the 2000-01 outburst of p 
Cas. The event was preceded by a strong collapse of 
the atmosphere which reversed to expansion and subse- 
quently cooled the entire atmosphere down by more than 
3000 K. The larger brightness decreases in R CrBs com- 
pared to the outbursts of p Cas are likely due to more 
vigorous dust formation around the former stars. The 
enhanced dust formation may be caused by the large car- 
bon abundance and hydrogen deficiency of R CrB atmo- 
spheres. Whether or not there is dust formation during 
the outbursts of p Cas, it appears reasonable to attribute 
the sudden eruptions to a common physical mechanism 
opera ting in the atmosph eres of cool hypergiants and R 
CrBs. lLobel et alJ (|2003f) showed that the recombination 
of partially ionized hydrogen in p Cas, and the recombi- 
nation of ionized helium in R CrBs can effectively drive 
a runaway explosion over the time scales corresponding 
to the observed brightness declines of both eruptive vari- 
ables. The eruptions are triggered by unusual thermal 
conditions in F-type stars during phases of very strong 
atmospheric compression c ausing a dynamicall y unstable 
atmosphere l|Lobe]|l200H) . lLobel et aT] (|2003) proposed 
that the stellar wind and oscillations can suddenly acti- 
vate the trigger, which we observe is tightly correlated 
with the occurrence of strong emission line spectra. The 
study of strong transient emission lines in different types 
of luminous pulsating cool stars therefore provides im- 
portant clues to the physical mechanisms that produce 
the large diversity of their observed light curves. 

R CrB stars both in and out of the sudden brightness 
declines reveal many spectral similarities with p Cas. 
Most noticeably the photospheric lines of R CrBs ar e 
very broad as in F-type hypergiants ijPandev et al.1 2004) . 
The lines sometimes split around phases of maximum 
brightness, which appears to correlate with the radial 
velocity amplitude. In RY Sgr, having a pulsation ampli- 
tude ±20 kms -1 , the split lines occur regularly, while in 
R CrB, having amplitude of only ±3 kms -1 , they are not 
obser ved every pulsation cycle, except in the strongest 
lines (Kamesw ara Rao fe LambertllT9 97). p Cas has an 
intermediate velocity amplitude of ±8 kms -1 , between 
these two values. The brightness declines are preceded 
by an increase of t he pulsation amplitude, obse rved for 
example i n R CrB ( Kames wara Rao et al.lll999|) . and in 
p Cas bv lLobel et all 1)2003). Possibly, it is therfore no 
surprise that the emission line spectrum also strength- 
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ens during variability phases preceding the outbursts of 
R CrBs. The so-called El group in R CrBs represents 
transient sharp emission lines that are only briefly ob- 
served at the very onset of the declines. The El group 
may therefore correspond to our variable type III emis- 
sion lines in p Cas. On the other hand, the so-called E2 
group of emission lines has spectral characteristics iden- 
tical to our type II lines. These emission lines are sharp 
cores inside low-excitation absorption lines, and are per- 
manently blue-shifted by 3-5 kins -1 relative to the star 
velocity. They are apparently of constant flux and appear 
most prominently in the middle of the fast brightness 
declines when the stellar continuum flux s ignificant ly de- 
creases llCottrell fc LamberlJ|1982t iGonzalez et al .lll998t 
IKameswara Rao et alJ 119991 ISkulian fc Cottrelll I2002TT 
Sharp emission line components above the local contin- 
uum level are observed during brightness minimum in 
R CrB stars in forbidden lines like [Ca II] and in reso- 
nance doublet lines of K I (0.77 pm) and Na D (0.59 
^m). The narrow emission components however do not 
occur in Mg n h & k line s (at 0.28Mm) or Ca II H fc 
K (0 39 pm) (IClavton et alJl99l lK ameswara Rao et aLl 
1999), which again compares to the observations of p 
Cas. Near-IR Na I or CO emission lines have never been 
reported in R CrB stars, possibly because this spectral 
region has not thoroughly been investigated. The optical 
C2 bands however are observed to appear in emission at 
the onset of the brightness declines. They either remain 
in emission or return to absorption around minimum 
brightness, although a clear correlation with the pulsa- 
tion phase or phases of brightness declines has yet to be 
estab l ished iGon zalez et alJl!998|: IKameswara Rao et al.l 
119991: ISkulian fc Cottrelll I200|k The two most signifi- 
cant spectroscopic differences between R CrB stars and 
p Cas are the presence in the spectrum of the former and 
the absence in the latter of: 1) a strong He 1 1.083 pm 
line with a P Cygni profile (although with notable ex- 
ceptions for warmer T e g > 6000 K stars as V854 Cen), 
indicative of warm chromospheric w inds in R CrB stars 
l|Clavt,on. Geballe. fc Bianchi 1200.^1 . and of 2) a broad 
(FWHM >200 kms -1 ) emission component observed 
during R CrB brightness declines in forbidden and per- 
mitted strong li nes (such as Ha, He, and me tal resonance 
lines), which IKameswara Rao et all l)1999l) suggest can 
result from an accretion disk around a putative white 
dwarf companion star (but note that a hot companion is 
not observed in p Cas). 

5.2.4. Line Splitting and Pulsations 

There is ample evidence from the literature that 
other long-period cool variables such as Cephcids, 
RV Tau stars, R CrBs, and Miras show absorp- 
tion line splitting and emission lines above the con- 
tinuum level around phases of maximum brightness, 
in atomic as well as molecular lines of CO and CN 
llSanfordll952HKovtyukh et al.l2003llAlvarez et all200H 
I Mozurkewich et al.l 11987} | Kameswara Rao fc Lambertl 
Il997t iNowotnv et al.ll2005bD . The line splitting phe- 
nomenon is also obse rved in hot RR Lyr variables 
ijChadid fc GilletHl9 96). Based on the above discussion 
we think it is merely a matter of coordinating coeval opti- 
cal and NIR observations to firmly establish that phase- 
dependent spectral changes in absorption line splitting 
and emission line flux are related spectroscopic phenom- 



ena, because spectral variability is clearly linked with 
changes of T e s in cool pulsating stars. These spectral 
phenomena are wide-spread among pulsating stars of var- 
ious types in the cool part of the HR-diagram. They 
should therefore be attributed to radiative transfer ef- 
fects that result from atmospheric pulsations instead of 
expanding circumstellar gas shells that are simply pos- 
tulated by default for detailed calculations of mass-loss 
rates. 

5.3. Emission by the Running Shock Wave? 

Considerable theoretical work has been published to 
model emission lines in cool pulsating stars. Among all 
cool stars with split absorption lines, the Mira variables 
are possibly best studied with recent dynamic models 
of the atmosphere that quite well match the observed 
line profiles. In most studies however the apparent ab- 
sorption components of split absorption lines are still in- 
terpreted to result from separate detached layers (op- 
tically thick 'shells') with different Doppler velocities. 
The double line cores a re simply assumed to have a 
purely kinematic origin lIBessel. Scholz. fc Woodl Il99ll 
IWinters et alJ HH INowotnv et a,U I2005alhfl. The al- 
ternative interpretation we propose with a central emis- 
sion reversal inside the absorption line core is compu- 
tationally much more challenging and has only very re- 
cently been explored with radiative transfer calculations. 
iRichter et alJ l|2003|) consider specific spectroscopic ex- 
amples with Fe II and [Fe 11] emission lines that occur 
in some Miras only after exceptionally large brightness 
maxima (or equivalently very strong atmospheric com- 
pression). They demonstrate that these emission lines 
can be excited in a confined post-shock layer close to 
the stellar photosphere at 1-2 The post-shock gas is 
heated by the passage of a shock wave produced in the 
photosphere, which reaches this emission layer in half a 
pulsation cycle. For the shock to be sufficiently strong 
to produce emission lines, a minimum shock velocity of 
~20 kms" 1 is required in Miras. In an earlier study 
iWoitke. Goeres. fc Sedlmavrl l)1996|) already demonstrate 
that adiabatic cooling of shocked gas can yield suffi- 
ciently small gas temperatures for the nucleation of dust 
grains within the same close distance of 1.5-3 i?* from 
photosphere of R CrB stars. The local gas density and 
shock velocity are the prime factors that determine the 
location of the emission zone. These models support the 
small distances of only <1 i?„ we compute for the NIR 
CO emission line formation region in p Cas with a sim- 
plified atmosphere model in which the outer gas temper- 
ature structure is modified by a shock wave. 

Some difficulties however remain when trying to di- 
rectly apply the above models to the warm yellow hy- 
pergiants p Cas and HR 8752. It should be remembered 
that radiation pressure on dust around the dust conden- 
sation radius (at a few i?*) is considered to be the driving 
mechanism of winds and mass loss in all modern models 
of R CrB and Mira variables. In p Cas however, there is 
no evidence so far of significant dust nucleation close to 
the photosphere. The weak 9.8 pm silicate dust emission 
feature and the IR spectral energy distribution place the 
dust envelope (possibly resulting fr om the outburst of 
1945- 46) at 30-130 R« from the star l|.Iura fc KleinmanrJ 
1990). New IR photometry and spectroscopy after the 
2000-01 outburst would help to further determine the 
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properties of circumstellar dust in p Cas. 

Furthermore, some constraints from our NIR CO emis- 
sion models can be used to estimate the velocity of the 
shock wave near p Cas. When index 1 denotes the pre- 
shock (upstream) gas conditions T\, pi, and pi, and 
index 2 the post-shock (downstream) conditions (with 
jump conditions across a compression shock surface: 
T2 > Ti, P2 > Pi, and p-x > pi), the speed of the shock 
Ui can be expressed through the local adiabatic speed of 
sound and the Mach number: u\ — c\ - Mi. For a specific 
heat ratio of ideal monatomic gas of 7 = 5/3 (assuming 
hydrogen gas) to 7/5 (for CO gas with rotational de- 
grees of freedom) and the mean gas-pressure and -density 
structures of our atmosphere model f *l4.2p . we compute 

ci=a/7 ■ pi/pi = 7.5-7 kms -1 . The pre-shock Mach 
number Mi can be computed with the classical Rankinc- 
Hugoniot jump relations in polytropic gas which are a 
function of 7 and T 2 /Ti only. For T 2 /Ti = 4200/2500 
K (Fig. 0) we compute Mi = 1.30-1.36. We thus ob- 
tain a shock velocity of u\ ^10 kms -1 . This value is 
smaller than the 20 kms -1 threshold velocity required 
in the shock models to excite atomic emission lines, and 
yet they are observed every pulsation cycle in p Cas. 
10 kms" 1 is in fact only a lower limit, since T2 (the gas 
temperature of the post-shock region) can be much larger 
than 4200 K we adopt from the unperturbed model at- 
mosphere. On the other hand, the non-detection of dust 
formation in the vicinity of the photosphere of p Cas ap- 
pears consistent with the small velocity (or Mach num- 
b er) of the shock wave (s). 

I Alvarez et alJ 1)20011) conducted a large survey of the 
absorption line doubling phenomenon in long-period 
variables. They demonstrate that it is more frequently 
observed among stars with compact rather than extended 
atmospheres, and interpret it to result from stronger 
shock waves in the former stars. Hypergiants such as 
p Cas and HR 8752 with their enormous stellar radii of 
several hundreds of R© do clearly not follow this trend. 
On the contrary, HR 8752 has become more compact as 
T e s has dramatically increased over the past 30 years, 
but the split absorption lines are presently no longer ob- 
served in its optical spectrum compared to the 1970s. If 
the line broadening of these split lines were solely due 
to a superposition of two absorption components emerg- 
ing from two separate gas layers above and below the 
shock surface, predicted with models of Miras, the width 
of these lines would have to decrease with the disappear- 
ance of the shock waves in the atmosphere of HR 8752. 
This is however not observed because the photospheric 
lines typical of a hypergiant spectrum have stayed very 
broad. 

Our model for the central emission excited by a prop- 
agating shock wave possesses some difficulties as well. A 
shock wave propagating supersonically in the stellar rest 
frame should produce red-shifted CO line emission (pro- 
vided the emission is sufficiently optically thick), which 
we do not observe. Another problem is that a transient 
and spatially rather narrow atmospheric structure re- 
quired for the CO emission cannot readily account for 
the permanent excitation of other static (and apparently 
centered around stellar rest) atomic emission lines in our 
high resolution spectra. 

In summary, none of the available theoretical models of 



the atmospheric dynamics can coherently explain the ab- 
sorption line splitting and line emission phenomenon in 
pulsating cool supergiants. The field of research is mak- 
ing rapid progress however. The current investigation 
shows that the correct answer may be obtained with the 
development of a new generation of atmospheric models 
that couple the dynamics of the atmosphere with non- 
LTE excitation of spectral lines in shock waves. 

6. SUMMARY AND CONCLUSIONS 

We describe a high-resolution spectroscopic study of 
one of the most luminous cool stars in the Galaxy, the 
bright yellow hypergiant p Cas. Our observations cover 
a complete pulsation cycle in 2003-04 over a large wave- 
length range from the optical to ~4 pm. Near-IR spec- 
troscopic studies of yellow supergiants are scarce, and we 
demonstrate that the emission lines they harbor serve as 
powerful diagnostics of the important atmospheric pro- 
cesses of mass loss and pulsation in massive stars. In 
particular, we identify several prominent emission lines 
in the L band spectrum observed near maximum bright- 
ness, and present if-band spectra obtained during the 
last grand outburst of the enigmatic hypergiant in 2000- 
01. 

We primarily investigate the first overtone band of CO 
at 2.3 pm because it occurs prominently in emission dur- 
ing certain variability phases of p Cas. The NIR CO 
band is believed to originate from expanding shells and 
to trigger the formation of dust close to t he photosphere 
of R CrB stars. We confirm the findings of lShefferl (1993) 
that the CO band in p Cas is variable and consists of ab- 
sorption and emission line components. However, based 
on detailed comparisons with atomic lines and correlat- 
ing all spectroscopic observations with the corresponding 
pulsation phases, we argue that the 'split' absorption line 
profiles are not caused by two separately ejected shells 
during rare outbursts of the hypergiant, but rather are a 
commonly observed combination of a static narrow cen- 
tral emission line superimposed on the core of a broad 
absorption line. The absorption line exhibits Doppler 
displacements due to the atmospheric pulsations, while 
the emission component is always centered around the 
stellar rest velocity. The flux of the central emission line 
varies with the pulsations of the photosphere and rises 
above the continuum level around maximum brightness 
during phases of fast atmospheric expansion. After the 
phase of minimum brightness the emission flux decreases 
and the lines fade. Detailed comparisons of NIR CO 
line profiles with low-excitation atomic lines reveal that 
they are intermediate between the type II and type III 
lines of p Cas (according our terminology). While CO 
absorption components are always present and vary sim- 
ilarly to permanently split atomic absorption lines (type 
II), the large variability of the superimposed CO emis- 
sion components suggests they are more compatible with 
type III emission lines that can appear above the level 
of the stellar continuum flux. In other words, the excita- 
tion mechanism of prominent atomic emission lines and 
CO emission we observe in p Cas can have a common 
physical origin. 

The strong correlation of the NIR CO emission flux we 
observe with the radial velocity curve from photospheric 
pulsations signals that the CO emission lines form in the 
vicinity of the photosphere. We model the CO emis- 
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sion band spectra by inserting a cool gas layer of T g 
^2500 K in the kinetic gas temperature structure of a 
plane-parallel hydrostatic Kurucz model of the photo- 
sphere with T cS =7250 K and log g =0.5. We fit the 
observed emission flux around the CO bandhead, which 
is least influenced by absorption, with radiative transfer 
calculations to constrain the optical depth of the cool CO 
layer. We determine an atmospheric column mass along 
the line of sight of io _31±01 g cm~ 2 . It translates to a 
geometric distance for the cool CO layer of only ~5% of 
i?* above the photosphere. 

We further compare the evolution of NIR emission 
lines in p Cas with the spectra of other late-type super- 
giants. We consider cool Miras, the yellow hypergiants 
HR 8752 and IRC+10420, some less-luminous post-AGB 
stars, and their population II analogs of R CrB stars with 
large visual brightness declines. We observe a number 
of important spectroscopic similarities between all these 
cool variable star types. The most remarkable are: 1) 
when several epochs of observations are available the NIR 
CO emission flux is always observed to vary; 2) the CO 
emission flux is largest around maximum brightness, also 
observed in prominent atomic emission lines; 3) promi- 
nent emission lines and split absorption lines are more 
frequently observed in stars with larger pulsation ampli- 
tudes or with larger radial velocity amplitudes; 4) the on- 
set of (sudden) strong brightness decreases in p Cas and 
R CrB stars results from a stronger than usual compres- 
sion of the photosphere followed by a stronger than usual 
expansion and cooling of the entire atmosphere with the 
appearance of the type III variable emission lines we pro- 
pose are analogous to type El emission lines observed in 
R CrB stars; 5) during the fast brightness declines of 
outbursts the initial emission line spectrum is replaced 
by another type of emission line, the so-called E2 lines 
in R CrB stars, and in p Cas the type II emission line 
cores of split absorption lines. We propose that the lat- 
ter types of emission lines are always present but become 
more obvious during fast brightness declines because the 
mean photospheric flux suddenly decreases. 

In summary, we identify two groups of narrow emission 
lines with FWHM < 50 km s _1 that coexist in the spectra 
of low-gravity pulsating stars. The first group consists 
of atomic lines that appear to be intrinsically invariable 
(types II and IV in our terminology). The second group 
consists of atomic and CO emission lines that are intrinsi- 
cally variable (type III). Both groups are low-excitation 
lines with xi ow of the order of few eV having Doppler 
velocities of ^5 kms -1 around the stellar rest velocity. 
The excitation of both groups appears to be related in 
one way or another to the atmospheric pulsations. The 
best interpretation we can propose based on our spec- 
troscopic observations is that the first group is excited 
in two st atic conic regions ^10-100 i?* above the pho- 
tosphere l)Lobel et al .1120 03). where the supersonic wind 
of p Cas collides with circumstellar material (and possi- 
bly replenished by recurrent outbursts) forming a perma- 
nently standing shock wave. The second group of emis- 
sion lines is excited in a short-lived circumstellar cool 
gas layer in the immediate vicinity of the photosphere 
less than 1 i?* away, that is excited by the passage of a 
pulsation-driven shock wave. 

The following list of observations would be helpful to 
further verify our interpretations: 



1. very high-resolution spectroscopic observations (R 
> 100,000) of type II and III emission lines in p 
Cas to search for possible differences in radial ve- 
locity and line profile shapes between both groups. 
The observations could confirm if the emission lines 
form in two spatially separate regions in the cir- 
cumstellar environment of the hypergiant; 

2. high-resolution spectroscopic observations beyond 
3 /im to identify and investigate the properties 
of other prominent emission lines expected in this 
wavelength region; 

3. monitoring of the fundamental band of CO at 4.6 
pm to determine how far out the photospheric pul- 
sations and shock waves propagate into the circum- 
stellar environment; 

4. new observations of the silicate 9.8 /jm dust emis- 
sion feature after the outburst of 2000-01. Near- 
and far-IR photometry can establish if p Cas pro- 
duces dust during the large brightness decreases of 
outbursts. It can further link its remarkable atmo- 
spheric dynamics to that of other eruptive variable 
stars; 

5. high-resolution spectroscopic monitoring of NIR 
CO and split absorption lines in other late-type 
stars. In particular, observations of semi-regular 
pulsating variables (e.g. RV Tau, R CrB stars) can 
test the relative importance of the pulsations versus 
T e ff for sustaining the cool gas layers that produce 
the peculiar emission line spectra. 
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TABLE 2 

Unindentified emission lines in the L 

BAND SPECTRUM OF p CAS. 



# 


A a 


FWHM 


EW 


(cm) 


H 


Om) 


1 


3.6584 (± 0.0015) 


0.0016 


0.00005 


2 


3.6991 (± 0.0015) 


0.0022 


0.00009 


3 


3.7456 (± 0.0015) 


0.0013 


0.00012 


4 


3.9115 (± 0.0016) 


0.0016 


0.00008 


5 


3.9885 (± 0.0016) 


0.0023 


0.00011 


6 


4.0521 (± 0.0016) 


0.0013 


0.00007 


7 


4.0609 (± 0.0016) 


0.0017 


0.00004 



a Rest wavelength in air - observed wave- 
length corrected for the Earth's motion 
and for the heliocentric radial velocity of 
p Cas (— 47±1 kms -1 ); uncertainties in 
brackets - width of the resolution element. 
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Fig. 1. — CorMASS R ~300 spectra of p Cas compared to other yellow hyper- and supergiant stars. The spectra have been normalized 
at 1.6 (im and shifted vertically for clarity. Clipped out are the noisy regions due to strong telluric absorption. Identified are the most 
prominent absorption and emission lines. The inset panel shows an expanded view of the K band spectrum. N otice that the NIR, CO band 
is only present in p Cas. The line identifications are based on high-resolution atlases of the Sun and Arcturus IHinkle et al.ll2003) . 



Fig. 2. — SpeX R ~2,000 continuum normalized emission spectrum of p Cas. Marked are some of the prominent emission lines previously 
reported in the literature, as well as some unidentified lines in the L band discovered in this work. The strong absorption lines beyond 1.5 
pm are mainly due to hydrogen. 



Fig. 3. — CGS4 spectra of p Cas observed on two consecutive nights in the middle of the deep brightness decrease of 2000-01, compared 
to other bright supergiants ordered according to spectral types, and to the emission spectrum of 2004. The comparison spectra are offered 
in Wallace & Hinklc (1997). All spectra have been convolved with the same instrumental resolution of i?~3,200 (except for the top SpeX 
spectrum which is of i?~2500), continuum normalized, and shifted vertically for detailed comparisons of line depths. The left-hand short 
horizontal dashes mark the zero-flux level for each spectrum. The spectrum of p Cas in outburst best fits the early K-type supergiant 63 
Cyg, although with noticeable exceptions (see text). The two outburst spectra are shown to demonstrate that the discussed features are 
real and not due to noise. 



Fig. 4. — Crosses: the visual light curve of p Cas from AAVSO binned in 10 d averages. The size of symbols is propo rtional to the number 
of averaged datapoints. Dashed line: radial velocity curve from photospheric absorption lines of [Lobel et al.I fcOQjf) . The star velocity is 
— 47±1 kms -1 . Arrows show dates of Shcffer (1993) CO observations presented with the same color scheme on Fig. 5c. Upper panel: our 
low-resolution K band spectra and one high-resolution spectrum from Hinkle et al. (1981) of p Cas. Near-IR CO is observed in emission 
on the descending branch of the light curve when star begins to expand out of minimum radius. 



Fig. 5. — a): high-resolution spectra of p Cas of the first overtone Ai/2-0 CO band observed in emission in 2004 September (solid line), 
and in absorption in 1979 February (dotted line; from lLambert. Hin kc & Hal| 119811^ . b): Separate rotational CO transitions of 2004 
September are overplotted in the velocity scale, together wit h Na I and Br7 lines observed on the same night, c): Averaged line profiles of 
the R14 to R17 CO transitions from 1979 to 1989 observed by Shcffer 1 1993) compared to our observations in 2004. The spectral resolutions 
are 9 and 8 kms - per pixel, respectively. The star velocity of —47 kms - is marked with the vertical dashed line. The spectra have been 
continuum normalized and shifted in the vertical direction. The pulsation phases are roughly determined from the light curve in Fig. [I] 
with 0=0.0 corresponding to maximum brightness, and 0= 0.5 to minimum brightness. 



Fig. 6. — The spectral energy distribution in the low-resolution CorMASS spectra can be used to constrain the photospheric parameters 
of p Cas, together with the temperature-sensitive hydrogen lines and the gravity-sensitive Ca II lines. The graph shows the best-fit model 
spectra from the NextGen-giant grid to the spectrum of 2004 October. The strong calcium triplet lines favor the lowest gravity value of 
0.0, for which models with only T cf f < 6400 K are availa ble. Near-IR CO is not expected to form in this temperature range confirming our 
calculations using Kurucz model atmospheres (see i|4.2l . 



Fig. 7. — Upper panel: the atmospheric gas temperature structures used to calculate the model spectra. The horizontal axis is in the 
column mass scale m. The outer layers of the atmosphere correspond to smaller values of log m. The marked numbers on the atmospheric 
temperature profiles show the corresponding geometric scale expressed as the height above the photosphere (where by definition the optical 
depth t =2/3 and T g = T e ff) in Kq (upper row), and in percent of the photospheric radius (adopting i?*= 450 Rq for p Cas; lower 
row). Lower panel: CGS4 high-resolution spectrum of p Cas of 2004 September (solid black line) plotted against synthetic model spectra 
calculated with temperature structures in the upper panel. The model spectra are blue-shifted by 7 kms -1 to match the velocity of the 
observed CO bandhead. Dot-dashed spectra are calculated using Kurucz atmosphere models, while the long-, short-dashed and dotted 
models include a cool gas layer in the upper atmosphere at log m = —3, —3.5, and —4, respectively. The model adapted with a temperature 
minimum at log m=— 3 best fits the observed emission flux around the CO bandhead. None of the adapted models yield significant NIR 
CO absorption. 



Fig. 8. — Examples of the main types of atomic spectral lines observed in p Cas, arranged in progression from pure absorption lines (of 
hydrogen) to pure emissions lines (of sodium). Compare with the line profiles of several Av =2-0 CO transitions. The oscillator strengths 
and excitation energies of the lower and upper energy levels for a given transition are shown in the lower right corner. Lines beyond 1 pm 
are from the CGS4 echelle spectrum of 2004 September, while l ines shortward of 1 ( im are from the ECHLR spectrum of 2004 October. The 
spectra of Arcturus, the Sun, and the telluric spectra are from IHinkle et all I2003fl . The NIR spectra are corrected for telluric absorption; 
weak telluric lines remain in the optical. The spectral resolution of all the NIR spectra and of the optical spectrum of p Cas is ~40,000, 
whereas for the optical spectra of the Sun and Arcturus R ~150,000, which in all cases sufficiently resolves the line profiles for direct 
comparisons. 



Fig. 9. — High-resolution optical Sofin spectra of p Cas during the same pulsation cycle as our NIR, CO observations. The upper panel 
shows variations in the high-excitation Fe I 5572.84 A absorption line which probes the variable photospheric radial velocity and changes 
of T e ff . The violet line wing forms in the upper atmosphere and probes opacity changes in the supersonic wind of the hypergiant. The 
core of the neutral iron line remains single throughout the entire pulsation cycle, and increases EW when T e g decreases. The lower panel 
shows three examples of lines that develop emission cores around maximum brightness (type III lines in our terminology). The flux of the 
emission lines is clearly intrinsically variable, independent of the changes in the underlying photospheric absorption spectrum. 
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Fig. 2 
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APPENDIX 

SPECTRAL LINE IDENTIFICATION LISTS 

Table lA"3l provides wavelengths of spectral lines marked in the figures of this paper. Sources of wavelengths we used 
are: 

• Atomic Line List v. 2.04, compiled by P. van Hoof 
www . pa . uky . edu/^peter/ atomic 

• 1995 Atomic Line Data by R. L. Kurucz and B. Bell, Kurucz CD-ROM No. 23. Cambridge, Mass.: Smithsonian 
Astrophysical Observatory 

cf a-www . harvard . edu/ amdata/ ampdata/kurucz23/ sekur . html 

• CO Line List by R. L. Kurucz 

kurucz . harvard . edu/LINELISTS/linesmol/ coxx . asc 
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TABLE A3 
List of observed spectral lines. 



A (A) 



Element 



Figure No. 



5572.842 

5269.537 

5371.489 

6191.557 

6108.116 

6562.80 

6643.630 

6645.064 

6663.441 

7323.89 

7771.944, 7774.166, 7775.387 

7912.866 

8047.617 

8203.6 

8498.03, 8542.09, 8662.14 

9229.014 

10036.653 

10049.373 

10114.015 

10123.86 

10124.940 

10327.311 

10396.76 

10455.449, 10456.757, 10459.406 

10683.08, 10685.34, 10691.24 

10914.887 

10938.095 

12818.08 

14584.2 

15024.993 

15700.662 

15880.541 

16109.313 

16407.192 

16806.520 

17108.663 

17362.108 

20629.757 

21163.76 

21354.226 

21655.29 

21782.7 

21874.15 

21897.3 

22004.3 

22024.3 

22056.4 

22083.7 

22210.6 

22232.9 

22257.108 

22266.7 

22274.0 

22624.97, 22626.73 

22651.18, 22653.58, 22655.35 

22928.976 

22972.326 

23106.100 

23214.513 

23220.222 

23518.163 

23822.970 



Fe I 
Fe I 
Fe I 
Fe I 
Ni i] 
Ha 
Ni i] 
Eu II 
Fe I 
[Ca II] 
O I 
Fe i] 
Fe i] 

Paiim 

Ca II 
PaC 
Sr II 
Pa<5 
Fe I 
C I 

Si i 

Sr II 

Ti I 
S i 

C i 

Sr II 

Pa7 

Pa/3 
Bri im 

Mg I 

Brig 

Bri4 

Bri 3 

Bri2 

Brn 

Mg I 

Brio 

Fe I 

Al i 

Si i 

Br7 

Ti i 

Si I 

Ti I 

Ti I 

Sc i] 

Na I 

Na I 

Ti I 

Ti I 

Fe I 

Sc I 

Ti I 

Ca I 

Ca I 

12 CO Au 2-0 bandhead 

AJ 36-37 

AJ 21-22 

AJ 13-14 
12 CO Av 3-1 bandhead 
12 CO Av 4-2 bandhead 
12 CO Ay 5-3 bandhead 



8, 9 



1, 3, 4, 5, 6, 



2, 3, 
2, 3, 



1, 2, 3, 4, 5, 



a Rest wavelength in air. 



